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ABSTRACT 

Context. Stellar bars play an essential role in the secular evolution of disk galaxies because they are responsible for the redistribution 
of matter and angular momentum. Dynamical models predict that bars become stronger and longer in time, while their rotation speed 
slows down. 

Aims. We use the Spitzer Survey of Stellar Structure in Galaxies (S^'G) 3.6 pm imaging to study the properties (length and strength) 
and fraction of bars at z = 0 over a wide range of galaxy masses (M, 10* - 10 ''Mq) and Hubble types (-3 <T< 10). 

Methods. We calculated gravitational forces from the 3.6 pm images for galaxies with a disk inclination lower than 65°. We used 
the maximum of the tangential-to-radial force ratio in the bar region (Qb) as a measure of the bar-induced perturbation strength for 
a sample of ~ 600 barred galaxies. We also used the maximum of the normalized m=2 Fourier density amplitude (A““) and the bar 
isophotal ellipticity (e) to characterize the bar. Bar sizes were estimated i) visually, ii) from ellipse fitting, iii) from the radii of the 
strongest torque, and iv) from the radii of the largest m=2 Fourier amplitude in the bar region. By combining our force calculations 
with the Hi kinematics from the literature, we estimated the ratio of the halo-to-stellar mass (Mb/M,) within the optical disk and by 
further using the universal rotation curve models, we obtained a first-order model of the rotation curve decomposition of 1128 disk 
galaxies. 

Results. We probe possible sources of uncertainty in our Qt, measurements: the assumed scale height and its radial variation, the 
influence of the spiral arms torques, the effect of non-stellar emission in the bar region, and the dilution of the bar forces by the dark 
matter halo (our models imply that only ~ 10% of the disks in our sample are maximal). We find that for early- and intermediate-type 
disks (-3 < T < 5), the relatively modest influence of the dark matter halo leads to a systematic reduction of the mean 2b by about 
10-15%, which is of the same order as the uncertainty associated with estimating the vertical scale height. The halo correction on 2b 
becomes important for later types, implying a reduction of ~ 20 - 25% for T = 7-10. Whether the halo correction is included or not, 
the mean 2b shows an increasing trend with T. However, the mean decreases for lower mass late-type systems. These opposing 
trends are most likely related to the reduced force dilution by bulges when moving towards later type galaxies. Nevertheless, when 
treated separately, both the early- and late-type disk galaxies show a strong positive correlation between 2b and A™’'. For spirals the 
mean s x 0.5 is nearly independent of T, but it drops among SOs (» 0.2). The 2b and s show a relatively tight dependence, with 
only a slight difference between early and late disks. For spirals, all our bar strength indicators correlate with the bar length (scaled to 
isophotal size). Late-type bars are longer than previously found in the literature. The bar fraction shows a double-humped distribution 
in the Hubble sequence (~ 75% for Sab galaxies), with a local minimum at T = 4 (~ 40%), and it drops for M, < 10® If we 

use bar identification methods based on Fourier decomposition or ellipse fitting instead of the morphological classification, the bar 
fraction decreases by ~ 30 - 50% for late-type systems with T > 5 and correlates with My,IM,. Our ratios agree well with 

studies based on weak lensing analysis, abundance matching, and halo occupation distribution methods, under the assumption that 
the halo inside the optical disk contributes roughly a constant fraction of the total halo mass (~ 4%). 

Conclusions. We find possible evidence for the growth of bars within a Hubble time, as (1) bars in early-type galaxies show larger 
density amplitudes and disk-relative sizes than their intermediate-type counterparts, and (2) long bars are typically strong. We also 
observe two clearly distinct types of bars, between early- and intermediate-type galaxies (T < 5) on one side, and the late-type systems 
on the other, based on the differences in the bar properties. Most likely this distinction is connected to the higher halo-to-stellar ratio 
that we observe in later types, which affects the disk stability properties. 

Key words, galaxies: barred - galaxies: evolution - galaxies: structure - galaxies: dark matter - galaxies: statistics 


1. Introduction 

The lambda cold dark matter (ACDM) model predicts that galax¬ 
ies originate in dark matter haloes. In the first stage of their evo¬ 
lution, galaxies undergo a process of continuous mergers. As the 

* Full Tables |A. 1[ |A.2| and |A.3| the tabulated radial force profiles 
and the rotation curve decomposition model of each individual galaxy, 
are only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A-l-A/ 


redshift decreases, galaxy mergers become less frequent and the 
evolution becomes internally driven (the so-called secular evo¬ 
lution). 


Nearly two thirds of all galaxies in the nearby Universe have 
bars (e.g. de Vaucouleurs et al.|1991[|Kna pen et al.|2000[|Why te 
et al.|2002[|Laurikainen et al.|2004a[|Menendez-Delmestre et al. 


2007|l. Approxim ately one half of these are strongly barred (e.g. 
de Vaucouleurs||T963[ Sellwood & Wilkinson||199'^ . Although 
early studies suggested a nearly constant bar fraction at least up 
to z = 1.2 (e.g.|Elmegreen et al.|2004|l, it has become clear that 
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bar fraction increases with time and is around 20% at z = 0.84 
( |Sheth et al.|2008||Nair & Abraham|2010| l. From all these studies 
it is known that the bar fraction depends on colour, stellar mass, 
and bulge prominence. 

Bars continuously interact dynamically with the other struc¬ 
ture components of galaxies, such as the underlying disks, the 
bulges, or the dark matter halos. This causes them to evolve 
over time (e.g. Hernquist & Weinberg|1992| Debattista & Sell 


morphology; |Danby||1965[ |Laurikainen et al.]|2007[ [Martinez^ 
Valpuesta et al.|2007|. Bar strengths, that is, the ratio of the tan¬ 


gential force to the mean axisymmetric radial force field, can 
be computed from the gravitational potentials inferred from op¬ 
tical and infrared images (Combes & Sanders 198 1[ Buta & 


wood|2000 Athanassoula||2003f Korinendy & Kennicutt||2004t 


Block|200Tl ILaurikainen & Salo|2002[ |Laurikainen et al.|2002| 


Berentzen et al.|2006[ Villa-Vargas et al.|2009| l. Galaxies are not 
isolated entities in the Universe, they are for instance actively 
interacting with the surrounding extragalactic gas and with the 
small satellite galaxies. Bars play an important role in the secular 
evolution of disk galaxies ( |Athanassoula|20l3] l. Bars are known 
to be robust so that their effects cover a wide range of redshifts 


2004a y. There are alternative ways of estimating the strength 

of the bars photometrically, namely by calculating their intrin¬ 
sic ellipticity ( Martin|[l995[ Laurikainen et al.||2002|l and from 
the normalized Fourier intensity amplitude ( [Laurikainen et'af] 
2004b 2005| l. Another indirect method is provided by studying 
dust lanes, whose curvature is inversely proportional to the bar 
strength according to the studies of|Athanassoula|(|1992 1 (probed 


(e.g. [Shen & Sellwood|2004[ Villa-Vargas et al.|2010||Athanas- 

soula et al. 

2013[l. This is supported by studies showing that bars 

are typical 

y composed of old stars ([Gadotti & de Souza[[2006[ 

Sanchez-B 

azquez et al. 120 11b. 


in Knapen et al. 2002[ Comeron et al. 2009 [ Sanchez-Menguiano| 


et al.|2015|l, or by 2D decomposition, fitting generalized ellipses 


to the bar component and using the boxyness as a proxy of the 
bar strength (|Gadotti 201 l|l. Bar lengths have also been mea¬ 


sured by different methods in the literature (for a detailed review 


Bars are responsible for the redistribution of the angular mo 
mentum of the baryonic and dark matter components of disk 
galaxies (e.g. Weinberg|1985 |Athanassoula & Misiriotis|2002 


|Marinova & Jogee 2007|l. Bar^induced secular evolution is mani¬ 


see |Erwin||2005 

hereafter E2005) such as ellipse fitting to bar 

isophotes (e.g. V 

^ozniak et al.|1995 1, Fourier analysis of the im- 

ages based on the bar-interbar luminosity contrast or phase (e.g. 

Ohta et al.|1990 

Aguerri et al.|2000( [Quillen et al.[[1994| [Lau- 

rikainen & Salo 

[2002|l, and they have been estimated visually 


fested in many features of the disk, such as resonance rings (e.g. 


Schwarzl 19811 |Buta|| 19861 IButa & CombesH 19961 IRautiainen 

From N-body simulations it has been shown that bars form 

& Salo|2000l, or long-lasting spiral density waves (e.g. Toomre 

spontaneously in galactic disks (e.g. Miller et al.[[l970 Hohl 

1969[ Kormendy & Norman|1979 1. 

T97T1 Ostriker & Peebles|19731[Sellwood|1980||19811[Athanas- 

Bars also trigger gas inflow towards the central regions of 

soula & Sellwood|1986 1, provided that the disk velocity disper- 


the galaxies, which leads to nuclear starbursts. They drive the 
secular evolution of bulges and cause the kinematic heating of 
the inner disk (e.g. [Combes et al.|1990| l. The causality between 
the bar-driven gas inflow and the fuelling of the active galactic 
nucleus (AGN) (e.g. [Shlosman et al.|[l989|l has been a much 


sion is not too high and the contribution of the disk to the to¬ 
tal force field is considerable. A key process in the bar-induced 
secular evolution is the interaction between the bar and the dark 
matter halo, which is expected to take place in resonance regions 
that are sufficiently filled with the baryonic matter ( Athanas^ 


debated topic during the past decades (for a review, see Jogee[ 

soula||2003|l. Theoretical models (e.g. |Villa- Vargas et al.| 20101 

2006[l. Laurikainen et al.[( 2004a i observed a higher fraction of 

Athanassoula et al.||2013|l predict that when the bar evolves, it 


bars in Seyfert galaxies than in the non-active counterparts; this 
consistent with other studies at near-infrared (e.g. [Knapen 


et al. 


et al. 


2000 


Laine et al.|2002 and optical wavelengths (e.g. Hao 


20091. [Lee et ah ([2012 1 recently claimed that the connec 


becomes longer, stronger, and slower. The last of these predic¬ 
tions is not always confirmed with the bar pattern speed obser¬ 
vations (e.g. Rautiainen et al.|2005| 2008|l, or at most, contradic¬ 


tion between the bar fraction and the nuclear activity only de¬ 
pends on the fact that AGN-host galaxies are on average redder 
and more massive, since no correlation is found when the colour 
or the stellar mass of the host galaxy are fixed. Likewise, [Cister-[ 
nas et al. (|2013[l did not observed any correlation between the bar 


tory observations are shown. In fact, direct Tremaine-Weinberg 
measurements (Tremaine & Weinberg[|1984]l indicate fast bars 


(jAguerri et al. 2bl5|l, but in this study our aim is not to solve this 


strength and the degree of nuclear activity as determined from 
Chandra X-ray observations, which suggests that other mecha¬ 
nisms are required to fuel the AGN (see also the discussion in 
[Hao et al.|2009[ and references therein). 

Bars themselves evolve, which among other things can man¬ 
ifest as buckling in the vertical direction that results in ver- 


Combes et al.[ 1990[|Raha et al.j 

19911 lAthanassoula & Misiri- 

Otis [2002 

[Debattista et al.[|2004 

Martinez-Valpuesta & Shlos- 


( [2007| l, [Laurikainen et al.| ( [2014[ l and [Athanassoula et al. ( 2014 1 


proposed that the so-called barlenses ( Laurikainen et al(] 2011 
are boxy/peanut bulges seen face-on. For a review of the proper- 


problem. Instead, we estimate bar strengths or the prominence 
of bars using different methods for a representative sample of 
galaxies in the nearby Universe that cover a wide range of galaxy 
masses. In a forthcoming paper these measurements will be con¬ 
nected to bar pattern speed estimates for the same galaxies. 

Infrared (IR) observations, where the dust absorption is 
lower than in the visual wavelengths, are good tracers of old 
stellar populations that probe the underlying mass distribution in 
galaxies ( [Eskridge et al.|2000] l, which makes this spectral range 
well suited for the study of stellar structures. For this reason, 
the mid-IR im aging of the Spitze r Survey of Stellar Structure in 


Galaxies (S^G, Sheth et al. 2010 1 , which includes 2352 galaxies 


ties of boxy/peanut bulges, the reader is refereed to Laurikainen 
[& Salo[ ( [2015[ ) and [Athanassoula|P015| l. 

Many studies have been devoted to characterizing bars from 
an observational point of view. Bars can have either flat or expo¬ 
nential radial surface brightness profiles with respect to the sur¬ 
rounding disk ( Kormendy|1982[ Elmegreen & Elmegreen 1985| l, 
and a few of them present blobs at the end of the bar (ansae 


in the nearby Universe of different masses and morphological 
stages (T) from which we take our sample, is ideal for studying 
the properties of stellar bars. 

In addition to ellipse fitting, we use Eourier decomposition 
as a method for estimating bar strengths. We evaluated the grav¬ 
itational potential from the flux in the 3.6 pm images and cal¬ 
culated tangential forces normalized to the axisymmetric force 
field. In these calculations we took into account what is learned 


from the previous studies using a similar method (e.g. Buta & 
Block 2001 [[Laurikainen & Salo|20021 [Laurikainen et al.|2002| 


2004a 


b|[20b51[Salo et al.|2010[): we included the effects of the 
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bulge while de-prqjecting the images and considered the effects 
of spiral arms. We also estimated the effect of the dark matter 
halo and the possible effects of the vertical thickening of the in¬ 
ner part of the bar associated with a boxy/peanut bulge. We re¬ 
peated the calculations with mass maps in which the non-stellar 
contaminants are eliminated, taken from Querejeta et al. (2015 i. 

The paper is organized as follows: in Sect. 2, we present the 
data used in this study and the sample selection. In Sect. 3, we 
provide a thorough description of the way in which we calculated 
bar strengths and the stellar contributions to the circular velocity, 
as well as the different tests that we carried out to probe possible 
sources of uncertainty in our methodology. In Sect. 4, we explain 
the different bar size measurements. In Sects. 5 and 6 we analyse 
the bar fraction and the various bar length and strength param¬ 
eters as a function of Hubble stage and galaxy mass. In Sect. 7 
we discuss the evidence provided by the different measurements 
for the growth of bars over a Hubble time. Finally, in Sect. 8 we 
summarize the main results and the implications for the secular 
evolution of disk galaxies. 


2. Data and sample selection 

2.1. S^G sample 


The 2352 galaxies of the S^G survey ( |Sheth et^ 2010 1 were ob¬ 
served in the 3.6 pm and 4.5 pm filters with the InfraRed Array 
Camera (IRAC; [Fazio et al.|2004|l, ins talled on-board the Spitzer 
Space Telescope ( Werner et al.|2004) l. The sample was defined 
based on data gathered from HyperLEDA ^Paturel et al. 2003 i 
with the following criteria: 


- Radial velocity Vradio < 3000 km s ', which corresponds to a 
distance of < 40 Mpc. 

- Galactic latitude \b\ > 30°. 

- Total blue magnitude corrected for internal extinction, 

<15.5 mag. 

- Blue light isophotal angular diameter D 25 >1'. 


All Hubble types (T) are included, although the requirement 
of HyperLEDA Viadio measurement means that many gas-poor 
early-type galaxies are missing from the original S^'G-sample 
(the observatio ns are currently e xtended to include the missing 
early types, see Sheth et al. 2013 1 . The S"'G sample covers a wide 


range of stellar masses that span over several orders of magni¬ 
tude. 

The surface brightness depth is pj,,(,^^(AB){1(t) ~21 mag 
arcsec'^. As the mid-infrared images are barely affected by dust 
absorption or polluted by star formation (little contamination 
from hot gas and PAHs; |Meidt et al.| 2012[ |Zibetti & Groves] 
2011|l, they trace the old stellar population in galaxies well. Pos¬ 


sible contaminants are controlled in the mass maps in which they 
are eliminated ([Querejeta et al.|2015|l. 


2.1.1. Selection of sub-samples 


The selection of barred galaxies is based on the classification 


in 

;4 


Buta et al. (2015 B2015 hereafter), made for the complete 


S^'G. We first selected disk galaxies (ellipticals in B2015 were 
excluded) with an inclination i < 65°, which leaves 1345 galax¬ 
ies, for which the gravitational potentials and the stellar com¬ 
ponent of the circular velocity were obtained. This inclination 


upper limit is very similar to that (i - 60°) used in Comeron 


I ' We acknowledge the usage of the database (http://leda.univ- 
lyonl.frl. 


jet al.j ( |2014j l and recommended in jZou et ak] P014| ). Of these 
non-highly inclined galaxies, 860 are barred according to B2015, 
forming our sample of barred galaxies for which bar lengths are 
measured visually. This means that we did not use any intrinsic 


bar ellipticity cut-off for the sample selection (as in e.g. Abra- 
jham et al.|1999| l, which means that oval-like structures are also 
included in our analysis as if they were normal bars. 

To calculate the bar force and determine bar lengths and 
shapes through ellipse fitting, we tried to cover the largest possi¬ 
ble number of barred galaxies within S^'G. However, some faint 
dwarf, irregular and/or gas-rich late-type galaxies have peculiar 
bars that are complicated to measure. Another reason to exclude 
some galaxies are offset bars, that is, bars whose centres are dis¬ 
placed with respect to the centre of the galaxy (NGC 1345 is 
a good example). These are typical among Magellanic and ir¬ 
regular galaxies (8 < T < 10). Of the 860 barred galaxies in 
B2015, 654 have bar sizes, axial ratios, and position angles es¬ 
timated from the isophotal profiles; bar force calculations based 
on Eourier decomposition and gravitational torques are made for 
599 systems. 


2.2. S^G pipelines 

Eive pipelines (Pl-5) are behind the S^'G data release and subse¬ 
quent scientific products: 


PI produced science-ready mosaics from the raw data. The 
pixel scale of the 3.6 pm mosaics after the image proc essing 
and reduction is 0.75"pixel ' (for further details see Sheth 
et al.|2010 1 . The EWHM of the images 


2 . 1 ' 


P2 created object masks automatically using SExtractor 
( jBertin & Arnouts|1996 1, which later went through a human- 
supervised verification and edition by hand (see also P4). 

P3 ( jMunoz-Mateos et al.||2015] describing also PI and P2) 
provides RC3-type parameters, such as the total magnitudes 
at the two S^'G bands and the isophotal radii at the sur¬ 
face brightness 25.5 mag arcsec^^ (R 25 . 5 ) obtained from the 


3.6pm images. Based on the calibrations from Eskew et al. 
P012| E2012 hereafter), total stellar masses were derived 
from the 3.6 and 4.5 pm fluxes in the following manner: 




( 1 ) 


where D is the distance to the galaxy in Mpc and F is the 
galaxy flux in MJy. 

- P4 ( jSalo et ar]|2015| ) is dedicated to the decomposition of 
the two-dimensional light distribution into different struc¬ 
ture components such as bulges, disks, bars, nuclear point 
sources, and various other disk components. The decompo¬ 
sitions were carried out using GALEIT ( jPeng et ^|2010| l 
with the help of IDE-based visualization procedures. In P4 
ellipse fitting was also made using the 3.6 pm images, and 
the disk orientation parameters were estimated based on the 
ellipticity and position angle of the outer isophotes (PAouter, 
foutei- = 1 - (^/a)outer, where a and b refer to the semi-major 
and semi-minor axes, respectively). The inclination of each 
galaxy was obtained under the assumption of an infinitesi¬ 
mally thin and intrinsically circular disk: i - cos“'(b/a)outer- 
The inclinations of the galaxies were visually checked by 
de-projecting the galaxy images and ensuring that no stel¬ 
lar structures appeared artificially stretched in face-on view. 
These P4 orientation parameters were used to de-project the 
galaxy images and also to convert ellipse fit parameters mea¬ 
sured in the sky plane to the intrinsic disk plane values. 
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Fig. 1. See next page for caption. 
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Fig. 1. Example of force calculation for the barred galaxy NGC4548. Panel A: The de-projected image of NGC4548 in magnitude scale with 
a range of 17-25 The blue dotted line indicates the visually estimated length and position angle of the bar. The light green ellipse 

corresponds to the bar radius from the isophotal fit (maximum ellipticity). The purple and light blue circles have radii equal to rA 2 and rqb, 
respectively. Panel B: Image of the galaxy after subtracting the axisymmetric m=0 component, together with contours tracing the force ratio 
Ft/ (Fr) (separated by 0.1 intervals). The dashed lines indicate the regions where the tangential forces change sign. The outer circle delimits a 
region of radius twice the size of the bar (the same as in Panel A). Panel C: The Ft / (Fr) force map (butterfly pattern) of the galaxy. The length and 
ellipticity of the bar are traced with black lines. The inner dotted circle corresponds to the bulge effective radii from P4 decompositions. Panel D: 
The thick solid curve indicates the normalized tangential force amplitude {Qj in Eq.[^ calculated from the force maps, using the nominal vertical 
scale height estimated from the radial scale length. The dashed lines correspond to the assumed upper (red) and lower bounds (blue) for the scale 
height. The thin green line corresponds to the bar-only force profile (see the text). Panel E: The normalized m = 2 Fourier density amplitude vs. 
radius. The dashed lines show higher order even Fourier amplitudes. Panel F: The phase of m=2 Fourier amplitude vs. radius. Panel G: The solid 
line indicates the circular velocity curve calculated from the image (dashed lines correspond to the same upper and lower bounds as in panel D). 
The horizontal dotted lines indicate the maxima of the observed H i velocity amplitude (green) and that of our calculated stellar component of the 
velocity curve. The solid orange line indicates the slope inferred from the linear term of the third degree polynomial fit of the inner part of the 
rotation curve (dashed orange line). The dotted vertical orange line delimits the region where this fit is made, which is taken between the galactic 
centre and the radius of the maximum rotation within one fourth of R 25 5 - The dark green vertical dotted line corresponds to a radius of 2.2 times 
the disk scale length (radius of the maximum velocity of an exponential disk;||Freeman|197'0), while the vertical dashed line traces the bar length. 


- P5 ( [Querejeta et al.|[2015[ ) separates the old stellar popula¬ 
tion from the contaminant non-stellar emission. Using the 
3.6 jum and 4.5 fim images (see |Meidt et al.|2012| l, a statis¬ 
tical independent component analysis (ICA) technique was 
used to decompose the galaxy emission into old stellar pop¬ 
ulation light and non-stellar emission associated either with 
hot dust, PAHs, or asymptotic giant and red super-giant stars. 


We used the 3.6 ;um science-ready PI images. Compared to the 
4.5 pm band, 3.6pm images have a greater depth that reaches 
stellar mass surface densities as low as one solar mass per square 
parsec ( |Sheth et al.|2010 1. They also have a smaller FWHM than 
4.5 pm images. 

Before the force calculations, the sky background was sub¬ 
tracted, and the images were cleaned by filling the masked re¬ 
gions (with foreground objects or image defects) with values 
obtained by linear interpolation from the adjacent good-image 
pixels. Distances and /? 25.5 values are from P3. The total stellar 
masses are also from P3, computed using Eq.[T] For the parame¬ 
ters of the bulge and the disk scale lengths (/ir) and orientations, 
we used the P4-values. 

When converting the 3.6 pm flux to mass in the force cal¬ 


culations, we used the formula given in Munoz-Mateos et al. 
( |2013| l to estimate the stellar mass from the raw 3.6 pm absolute 
magnitude (M 3 . 6 AB): 


log(M*/Mo) = -O. 4 M 3 . 6 AB + 2.13, 


( 2 ) 


not available, we derived M/ from the B-band absolute magni¬ 
tude (e.g. Persic et al.|1996| l: M/ = 1.087 ■ (Mg + 0.38). 

We identified the bars based on the morphological classifica¬ 
tions by B2015, which were made visually using the 3.6 pm im¬ 
ages and following the comprehensive de Vaucouleurs-revised 
Hubble-Sandage system ( |de Vaucouleurs| |1959[ |Buta et al.| 
|2007] l. This catalogue uses the following notation for the promi¬ 
nence of the bar; SB, SAB, SAB, and SAB. The underline 
notation means that the galaxy is more likely barred or non- 
barred, hosting a strong or a weak bar. Magellanic and irregular 
barred galaxies (I) are also included in our analysis (Hubble type 
T=10). The maximum ellipticities and visual measurements of 
bars used in this study are taken from the catalogue of structures 
by |Herrera-Endoqui et al. ( 2015[ hereafter HE2015). 

The Cosmic Flows project provides a complete database of 
HI line profiles obtained at Green Bank in the USA and at Parkes 
in Australia ( |Courtois et al.|[2()09| |2011| l, from which we ob¬ 
tained the maximum circular velocity of our galaxies in a similar 
manner as for instance in|Zaritsky et al.|([2014|l: 


W^S:./(2sini), 


(5) 


where is the line width and i is the P4 inclination. For 
the galaxies for which Cosmic Flows data were not available 
(~ 40 %), we obtained the gas velocity amplitude from Hyper- 
LEDA, and corrected it to correspond to P4 inclination. 


assuming a mass-to-light ratio MjL - - 0.53 (E2012). In 

practice, Eqs.[T]and|2 are very similar, given the narrow range of 
[3.6] - [4.5] colours ( Munoz-Mateos et al.|2013| l. 


2.3. Other sources of information 

Absolute magnitudes in the I and B-bands (needed in Sect. |3.6| l 
were calculated from the apparent magnitudes available in Hy- 
perLEDA, which were corrected for Galactic (AG; Schlegel et al. 
|1998| l and internal extinction (Al; Bottinelli et al. 1995^ and K- 
correction (AB/;|de Vaucouleurs et al.|19761l as follows: 


niB = -AG-AI- AK(T) ■ v/10000. 


nil — ni'^’^ + 0.44(mB 


„leda\ 


(3) 


(4) 


3. Gravitational potential and calculation of bar 
forcing 

3.1. Method for caicuiating bar-induced tangentiai forces 

We inferred the gravity potentials from the 3.6 pm images and 
constructed maps of the tangential forces (Ft) normalized to the 
azimuthally averaged radial force (Fr). The calculations were 
made with the NIRQB-code (jUaurikainen & Salo||2002[ here- 


after LS2002), which is based on the polar method of Salo et al. 
\\999) . Figure shows a typical example of force calculations 
in a barred galaxy, with a well-defined butterfly pattern ( |Buta| 
& Block|2001|l in the torque map, which is roughly symmetric 


with respect to the bar major axis (with some twists caused by 
the spiral arms). From the torque maps we constructed the radial 


profile of the normalized tangential force amplitude (Combes & 
|Sanders|1981| l: 


where v is the heliocentric radial velocity (see HyperLEDA doc¬ 
umentation for further details). When /-band photometry was 


max (IFxCr, 0)|) 

<|FR(r,0)|) 


(6) 
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The maximum of Qj in the bar region characterizes the bar 
strength, called in the following. The radial distance where 
the maximum occurs is denoted as rgb. 

Instead of computing the gravitational potential directly from 
the image pixels, as is done in the Cartesian method (Quillen 
et al. |1994| l, the de-projected image was first decomposed into 
polar coordinates ( Salo et al.|1999[[2010| LS2002). This is based 
on the azimuthal Fourier decomposition of the surface densities 
in different radial zones: 


2(r, <p) = /o(r) 


1 + ^ A„(r)cos[in((fi - (pm(f))] 


m=l 


(7) 


The amplitude of each Fourier component, A„„ was tabulated 
separately as a function of radius (we note that amplitudes are 
normalized with respect to the azimuthally averaged m - 0 
component). Altogether, the summation was made over the even 
Fourier modes m-Q-2Q. The main modes in the bar region are 
m=2, 4, 6, and 8 ( |Ohta|[T996l LS2002). The maximum of the 
m - 2 Fourier amplitude in the bar region was used as another 
proxy of the prominence of the bar (denoted by A™“). The radius 
where this maximum is found is called r\ 2 - 

The advantage of the polar method is that it effectively sup¬ 
presses possible spurious maxima that may arise in the noisy 
outer parts of the images ( |Salo et al.|2004 1 . The potential at the 
equatorial plane is obtained by 


cl)Jr,<^,z = 0) = -G 


poo p2n 

I r'dr' I : 

Jo Jo 


<fi')g(Ar)d<fi', (8) 


where Ar^ - r'~ + - 2rr'cos(0' - 0), G is the gravitational 

constant, and g(Ar) is a convolution function including the pre¬ 
tabulated integration over the vertical direction (see LS2002). 
We note that the azimuthal integration can be made with FFT for 
each m while integration over radius is done by direct summa¬ 
tion. The assumptions made in our nominal force calculations 
are that (1) the mass-to-light ratio (M/L) is constant, (2) the 
vertical profile follows an exponential law, (3) the scale height 
is constant over radius, (4) the scale height is tied to the disk 
scale length. The effect of varying the assumptions (1) and (3) 
are studied below (Sects. |3.3| and |3.5| l. The role of the func¬ 
tional form of the density distribution (2) was explored in de¬ 
tail in LS2002, where it was concluded that Qb is not sensitive 
to the detailed vertical profile as long as the same dispersion is 
assumed (for the exponential law V< > = 2h^). 


3.1.1. Bulge stretching correction 


When necessary, a bulge stretching correction was also applied 
in the force calculations to prevent obtaining artificially strong 
tangential forces caused by the bulge appearing elongated af¬ 
ter de-projection (as in e.g. Laurikainen et al.||2004b| l. Thus for 
a prominent bulge, its flux is subtracted from the image before 
de-projection of the image to face-on orientation. The obtained 
force then represents the disk-only contribution. The contribu¬ 
tion of the bulge force is added, calculated by assuming that the 
bulge has a spherically symmetric intrinsic light distribution (in 
which case it is unambiguously determined by its projected light 
distribution). To obtain the flux of the bulge, the model parame¬ 
ters were taken from P4-decompositions (the fitted bulge model 
also corrects for the smearing due to the seeing) in which the 
model assumed elliptical bulge isophotes. We used an equiva¬ 
lent spherical radius that kept the fitted bulge flux constant. The 
typical result of the bulge correction is to remove the artificial 
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Fig. 2. First and second rows: The radial force profile of NGC 3953 and 
NGC 3504 with and without spherical bulge correction. The former (red 
line) represents the Qt profile when the bulge flux is subtracted from 
the image before de-projection of the image to face-on orientation and 
added back afterwards for the potential calculation. For the first galaxy, 
the global maximum of which corresponds to Q\y, is slightly over¬ 
estimated if no correction is considered, but remains within the limits 
determined by the disk thickness uncertainty (black dashed lines). The 
visual estimate of the bar length and the effective radius of the bulge are 
indicated with vertical dashed lines. Third and fourth rows: Compari¬ 
son of (2b measurements with and without bulge stretching correction 
for the galaxies with bulge in P4, with the data binned based on the disk 
inclination and the bulge-to-total ratio (P4). The strength of the correc¬ 
tion is also studied in terms of the location of rQb with respect to the 
effective radius of the bulge from P4 {rfj, separating the galaxies ac¬ 
cording to their morphological class. The outliers in the third row are 
also indicated with larger symbols in row 4. 


inner peaks of the profile in the bulge region (see the exam¬ 
ples in the two upper frames of Fig. |^. Naturally, the assump¬ 
tion of spherical bulges does not work well for non-classical 
disky bulges: the force calculations were thus also re-ran with¬ 
out any bulge correction, that is, assuming that the bulge has the 
same flattening as the disk. A third possibility, that the bulge 
represents the inner boxy/peanut part of the bar, is briefly ad¬ 
dressed in Sect. 3.5 To avoid over-correcting the images, the 
P4-decompositions were visually inspected to exclude highly ex¬ 
tended bulge fits. 

Altogether, the bulge correction and its uncertainties are 
expected to affect only galaxies that have relatively extended 
bulges compared to the bar size ( Laurikainen et al.|2004b l. In¬ 
deed, as shown in Fig. ^ the strength of the bulge correction 
is not very relevant in the statistical sense. The few outliers in 
























































S. Diaz-Garcia et al.: Characterization of galactic bars from 3.6;rm S^'G imaging 


the plot correspond to galaxies with very small bars, with rQb/re 
close to unity, where is the effective radius of the bulge. In 
total, 251 of the galaxies in our sample of barred galaxies have a 
bulge fitted in P4, of which 145 galaxies present central peaks in 
the force profiles that are due to bulge de-projection. However, 
only in 65 cases did bulge stretching change the value no¬ 
ticeably. For these we obtain median of 0.9 + 0.02 

and a mean |AQbl of 0.03. Overall, the role of bulge correction is 
small, and the different treatments used are expected to delimit 
the uncertainty in bulge contribution to Qh quite well. 

3.1.2. Vertical scale height 

We estimated the vertical scale height based on the empirical re¬ 
lation from |de Grijs| ( |1998| l, in which it is shown that the ratio 
of the disk vertical thickness to the disk scale length depends on 
the Hubble type. The morphological stage binning and the corre¬ 
sponding mean values (and range) are the same as in|Laurikainen] 
|etar] ( |2004bl l: /jr/Zz^ = 4(1 -5)if T < 1, 5 (3 - 7) if T e [2,4] 
and 9 (5 - 12) if T > 5. To estimate the largest uncertainties on 
the forces, the gravitational potentials were also calculated using 
the highest and lowest values for each bin (the number listed in 
the ranges inside the parenthesis). 

When /zr values were not reliable (129 cases in the sam¬ 
ple of 1345 non-highly inclined disk galaxies have quality flags 
lower than 5 in P4), we assumed that the vertical scale height 
scales with the disk size as - O.lr^ , where is the K- 
band isophote radius at 20 mag arcsec^^ from 2MASS0 This 
empirical relation, also used in |Salo et al.| ( |2010| ), was found in 
[Speltincx et al.| ( |200^ to approximate the de Grijs relation in 
an adequate manner R Finally, we used a25/l-5 in substitution 
of when neither nor /zr values were available (17 galax¬ 
ies). Here a 25 denotes the radius of the Z?-band surface bright¬ 
ness isophote of 25 mag arcsec^^, taken from HyperLEDA. This 
relation between these isophotal B- and K-band radii is purely 
empirical. 

The uncertainties of our force calculation method are sum¬ 
marized by |Laurikainen et ah] ( |2004b| , where it was shown that 
the largest expected contribution (10-15 % in the value of 2b) 
comes from the uncertain vertical thickness. On the other hand, 
the functional form of the vertical density distribution (LS2002), 


superposition of spiral arms (|Buta et aL|2003 

1 , the position an- 

gle of the bar relative to the line of nodes ( 

Buta et al. 2004 


BLS2004 hereafter), possible contribution of the dark matter 
halo (BLS2004), or the radial variations in the vertical scale 
height (LS2002) were estimated to affect Qt, by only about 5%. 
However, since our current sample extends to galaxies of much 
later type, a more thorough estimation of the effects of the dark 
matter halos is needed (see Sect. |3.6| l. 

3.2. Circular velocity curves 

We used the mass-to-light ratio at 3.6 /zm by E2012 (T 3 = 

0.53) and assumed it to be constant throughout the galaxy to 

^ 2MASS is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Tech¬ 
nology, funded by the National Aeronautics and Space Administration 
and the National Science Foundation. 

^ The bar torque parameter (Qg) was also computed based on this re¬ 
lation ( jSpeltincx et al.|2008[ for all the barred galaxies in our sample 
with an identifiable maximum torque. It is not used in the discussion 
except for Figs.|^and[^ Both approaches for the disk thickness deter¬ 
mination give very similar Qh values. All the measurements are listed 
in Tablejg 


obtain the stellar contribution to the circular velocity curve: 


F3.6/jm(f) = <ER(r)) r. 


( 9 ) 


where r is the galactocentric radius and Er is the radial force cal¬ 
culated for MjL - 1. We calculated the maximum of these ro¬ 
tation curves (F™^^) and the radius where V’ 3 . 6 /jm peaks (r™j]jjj) 
for all the 1345 disk galaxies with inclinations lower than 65°, 
regardless of whether they host a stellar bar or not. Uncertain¬ 
ties in the determination of Vs e/im are partly influenced by the 
disk scale height uncertainty (see Fig.[^. However, unlike in the 
calculation of normalized tangential forces, the uncertainty in 
the mass-to-light ratio (~ 30% according to E2012) is dominant. 
Variatio ns in Y 3 gmn are as sociated with the [3.6] - [4.5] colour 
(E2012, Meidt et al.|2014 1, but the colour differences in the 
sample are small enough (average of-0.41 +0.08 mag (AB), ac¬ 
cording to |Munoz-Mateos et al.|2013] l, to safely apply the same 
mass-to-light ratio to all Hubble types without any systematic 
error. 

Furthermore, we obtained the inner stellar velocity gradient 
of V 3 . 6 /jm, dRF 3 . 6 ;jm( 0 ), by fitting the inner rotation curve with a 
polynomial function of order m - 2> and taking the linear term 
as an estimate of the inner slope (in a similar manner as in|Lelli| 
|etal.|2013» . 


3.3. Qh from 3.6 fj.m images and from P5 mass maps 

In spite of the azimuthal smoothing in our polar method and 
the use of near-IR imaging, non-stellar emission in H ii regions, 
which also appears along the bar, may affect the Qj radial pro¬ 
files. This is the case particularly for the low-luminosity late-type 
galaxies, in which this dumpiness induces local maxima that can 
be mixed with the bar-induced amplitude maxima. 

To test the impact of non-stellar contaminants in our esti¬ 
mates, the force calculations made for the direct 3.6 pm images 
were repeated for a random sample of 72 barred galaxies using 
the P5 mass maps ( Querejeta et al.|2015[ ), including only galax¬ 
ies with reliable ICA analysis. For 34 of these galaxies, it was 
possible to reliably identify Qh . Two examples of force profile 
comparisons are shown in Fig. P5 indicates a large contri¬ 
bution of non-stellar contaminants for NGC 1566, in particular 
associated with the spiral arms, whilst for NGC 0150, the role 
of contaminants in the bar strength calculation seems less im¬ 
portant. Figure also shows the comparison of bar strengths de¬ 
rived from the direct images and from the mass maps for our 
sub-sample. The deviation in this relation is ~ 10-15 % (see 
also Fig.|^. The median 2b(ICA)/Qb(3.6/zm) is 1.06 ± 0.03, 
and the mean |A2bl = l2b(ICA) - Qb(3.6/zm)| is 0.03, which is 
smaller than the estimated uncertainty related to the thickness of 
the disk. Although some individual galaxies can be affected by 
the non-stellar contaminants, there is no systematic difference in 
the two Qh values. 

It might naively be expected that using mass maps would 
increase the number of low-luminosity galaxies for which bar 
strengths can be reliably estimated. However, even though the 
ICA method eliminates non-stellar components that are associ¬ 
ated with star-forming regions, for most of the low-luminosity 
galaxies in our sub-sample the Qj profiles also remain noisy 
when mass maps are used. 

In Appendix [C] we also discuss the effect of non-stellar con¬ 
taminants on the Fourier amplitudes of our bars and on the cal¬ 
culation of the stellar component of the circular velocity curves. 


Article number, page 7 of 38 










































A&A proofs: manuscript no. astroph_2015_AA_diazgarcia_simon 


NGC 1566 


Disk plane (Inclination ~32*) 




Qt, (3.6 pm) 


Fig. 3. Panel A: Radial force profiles derived from 3.6pm imaging 
(solid black line) and from P5 mass maps (red solid line) for the barred 
galaxies NGC 1566 and NGC 0150. For NGC 1566, which is represen¬ 
tative of the galaxies with the largest difference in the force profiles 
due to non-stellar contaminants, the Qj derived with the P5 mass map 
slightly differs from the raw Qj, hut the deviation is lower than the 
uncertainty associated with the disk thickness, whereas for NGC 0150 
the Qf, is roughly the same. We also show the “bar-only” Qj profiles 
using 3.6yum images for these two galaxies. Panel B: For a sample of 
34 hatred galaxies, we compare the gravitational torque parameter as 
measured from 3.6 pm and IGA corrected maps. The error bars are de¬ 
termined hy the disk thickness uncertainty. No systematic deviation is 
observed. 


3.4. Contribution of spiral arms 

The Qt profile of NGC 1566 (see the upper plot in Fig.[^ is rep¬ 
resentative of a barred galaxy with prominent spiral arms, pre¬ 
senting a local maximum due to the bar perturbation at a radius 
of ~ 30 arcsec, while the global maximum at ~ 85 arcsec corre¬ 
sponds to the spiral arms. 

We wish to determine the contribution of the spiral arms 

( 2003| l developed a 


to our bar force measurements. Buta et al. 


Fourier-based method for filtering the spiral atm contribution to 
the bar strength measurements in near-infrared images. They as¬ 
sumed that the A 2 (r) Fourier amplitude profile due to the bar 
alone is symmetric with respect to its maximum, rA 2 , if the 
effect of the spiral arms were suppressed, and used this sym¬ 
metry assumption to eliminate the spiral contributions from A 2 
(and likewise from A 4 ). However, this method is fairly time- 
consuming to apply and it is unclear how well it applies for the 
galaxies presenting inner non-axisymmetric structures (e.g. dou¬ 
ble barred galaxies) and for the early-type galaxies that have an 
ansae morphology. Moreover, using measurements from the OS- 
UBSG sample, [Laurikainen et al. ( 2007| l showed that the trend 
of 2b in the Hubble sequence was not affected by whether this 
spiral arm correction was applied or not. 


Salo et al. (20101 found a statistically significant correlation 


between the local bar-only forcing and the local spiral density 
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Fig. 4. First panel: We compare different models of a disk scale hight 
profile for NGC 4548. The profiles assume a Gaussian function, with the 
spread A defined in terms of the har size. The black horizontal dashed 
line represents the nominal thickness to this particular galaxy, while the 
uppermost and lowermost lines correspond to the limits for the disk 
thickness as determined in |Speltincx et al.| ( |2008^ . Second and third 
panels: Radial force profiles and rotation curves inferred from the po¬ 
tential associated with the distance-dependent exponential scale heights 
defined above (solid lines) and under the assumption of radially con¬ 
stant disk thickness (dashed lines). 


amplitude. In this study the spiral contribution to forcing was 
eliminated simply by setting the m > 0 Fourier density ampli¬ 
tudes to zero beyond a certain radius (rcut) (the polar method 
facilitates this separation). We followed a similar method here 
and chose rcut to equal the bar length (see Fig. 1 in |Salo et ah] 
pOlO I and Fig.[^in this work for an illustration of force profiles 
with and without spiral contribution). 

We denote as the value of the bar torque parame¬ 

ter after excluding the spiral arm contribution to the force maps 
in the bar region. In the statistical sense, bar-only force mea¬ 
surements do not deviate significantly from the raw Qt, v alues 


(^gbar iQ^ — 0.89; CT = 0 . 12 ), in agreement with 


Lau¬ 


rikainen et al. (2007 1 . To test the sensitivity of ^ to the 

assumed rcut, we also computed the bar-only profiles after set¬ 
ting Tcut = L2rbar, obtaining a mean gl’” of 0.97. 


3.5. Radial profile of the vertical disk thickness 

Our standard force calculations assumed a constant scale height 
throughout the disk. As mentioned above, LS2002 showed that 
shallow gradients in vertical thickness do not have a significant 
effect on 2 b compared to that obtained by assuming a constant 
average h^,. 

Here we determine the change in 2t for a more strongly 
varying disk thickness, for example in the presence of inner 
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itself (for the model with the smallest A, the force is enhanced 
by ~ 15%). The rotation curves derived from the azimuthally 
averaged radial forces were also compared and led to a similar 
conclusion: the resulting curves show very similar values for the 
four models (they drop by less than 10% at the thick bar region 
and they closely approach the rotation curve associated with the 
thickest disk under the constant /zj approximation). 

Our selection of /z™“ and /z™" implies a difference of a fac¬ 
tor of two between the inner and outer bar thicknesses, which 
is taken as an upper bound. To verify whether this limit is re¬ 
alistic, we examined the vertical profile of the edge-on galaxy 
NGC 4565 (see Fig. [5]l, wh ich is known to host a prominent 
peanut-shaped bulge Parvis|[r986| l. For this particular case, we 
estimated the disk scale height from the line-of-sight vertical dis¬ 
persion: 
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Fig. 5. First panel: Image of the edge-on galaxy NGC 4565, which 
hosts a prominent peanut-shaped bar, in magnitude scale with range 
17-24 yU 3 6^n,(AS). The galaxy has been rotated so that the major axis 
appears horizontal in the image. The blue box roughly covers the bar re¬ 
gion and extends vertically as far as Zs, which is the height above which 
90% of the light emission comes from the thick disk component. Sec¬ 
ond panel: Disk scale height estimated from the line-of-sight vertical 
dispersion. The continuous, dashed, and dotted blue lines correspond to 
the resulting integrating up to Zs, 1.5zs . and 0.5zs, respectively. For 
the first case, the black dashed horizontal line shows the mean outside 
the range of the peanut-shaped bar. The horizontal lines trace the thin 
disk Izz estimated from the fits to the vertical luminosity profiles (red). 
The resulting Qj profiles from the different disk scale heights (solid 
blue, dashed black, and dotted red lines; with the meaning indicated in 
the legends) are shown in the lower right plot. 


thicker structures such as the boxy/peanut bulges. These are 
vertically thick inner regions of bars that have gone through a 
buckling episode (see simulation s of|Combes & Sanders [T981[ 
Combes et al. 1990[ Raha et al. |1991| l and are readily identi- 
fiable in observations of nearly edge-on galaxies Parvis|[r986[ 


|Lutticke et al.|2000a|); they appear as X-shapes in unsharp mask 
images ( Bureau et al.|2006} Laurikainen et al.|2014| l. Moreover, 
it has been suggested by [Laurikainen et al.| ( |2007[ ) that the round 
inner parts of bars seen in many more face-on galaxies, termed 
barlenses (Laurikainen et al. 2011 1 , might in fact be the same 
phenomenon. Recent support for this speculation has been given 
in [Laurikainen et Jr| ( |2014[ ) and jAthanassoula et al.| ( |20T4| l. 

To obtain a rough estimate of the possible effect of a thick¬ 
ened central bar component, we calculated the gravitational 
forces by assuming that the disk scale height declines radially 
from /z™“ = 0.15 ■ rKj,, to /z™" = 0.075 ■ rx^^ following a function 

/z,(r) = /z™" -H (/z”“ - /z™") ■ exp( - 1/2 ■ (r/A)^), (10) 

with the range A varying from 0.25rbar to rbar. In Fig.j^we assess 
the effect of the disk scale height variations on the measure¬ 
ment for the particular case of NGC 4548. For the studied range 
of As, we observe that the differences in Qo are of the same or¬ 
der as the uncertainty arising from the unknown disk thickness 


^ V< ^ . 


P S{x,z)z^dz 

‘7—Zs 

\ p S(x,z)dz 


( 11 ) 


where x and z are the axial and vertical distances, S is the flux 
density, and Zs is the height above which 90% of the light emis¬ 
sion comes from the thick disk (taken from |Comer6n et al.|2012] l. 
Varying Zs by a factor 2 has no significant effect on the final 
radial profile, as shown in the same figure. We obtain a mean 
/Zz value in the bar region (~ 7 arcsec) consistent with Com- 
jeron et ah] ([2012 1, who performed fits to the vertical luminosity 


profiles in the disk region outside the peanut-shaped bar area. 
More importantly, we observe a difference between the maxi¬ 
mum /zz at the peanut-shaped bulge and that of the surrounding 
thin bar/disk (< 50%) that is much lower than the above tested 
values. The same most likely holds for the true as a function 
of the galactocentric distance. 

In the same figure we also show the force profile of 
NGC 4548 assuming it has the same scale height as NGC 4565, 
under the assumption that the peanut-shaped part covers one half 
of the bar size ( jLiitticke et al.|[2O0Ob| l. These two galaxies sup¬ 
posedly have similar bar lengths, and both are classified as Sa in 
B2015. The resulting Qt barely deviates from the force profile 
under the assumption of constant thickness. 

Using an image taken from an N-body-i-SPH simulation, 
Fragkoudi et al. (2015|) concluded that Qo can be strongly over¬ 


estimated by not taking into account the boxy/peanut bulge ge¬ 
ometry (regardless of measurement uncertainties). Determining 
the scale height in the peanut-shaped region from observed non- 
highly inclined galaxy images is complex (as discussed also in 
this work). We have analysed different models of a thick in¬ 
ner bar (twice the thickness of the surrounding thin disk), with 
the functional form of the scale height radial profiles defined in 
terms of the bar length. It is probably safe to conclude that the 
uncertainty on Q\, coming from the assumption of a constant /zz 
is not larger than the uncertainty arising from the disk thickness 
determination itself. In addition, for a complete view of the pos¬ 
sible effect of the boxy/peanut bulge, projection effects on the 
analysis also need to be included (in a similar manner as in jSaloj 
jet al.|20()4] l, which can be particularly strong in the inner regions 
(as discussed in Sect. 3.1.1 without taking into account the de¬ 
tailed bulge geometry). 

3.6. Dark matter halos 

Our goal is to estimate the amount of dark matter halo in the 
disk regions of our sample galaxies and to quantify the dark halo 
contribution to the mean axisymmetric radial force field, and 
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Fig. 6. Upper frame: The halo-to-stellar mass ratio within the optical 
radius for the non-highly inclined spiral galaxies (r < 65°) in the S°^G 
sample with usahle Hi maximum velocities (1227 systems). Galaxies 
with > 400 km s“* (gas-poor galaxies with unreliable velocity 
measurements) are excluded from our statistics. Both barred and non- 
barred systems form part of this sample. Moderately inclined galaxies 
(30° < i < 50°) are plotted with orange filled circles, while more in¬ 
clined and face-on galaxies are displayed with orange X symbols: for 
both sub-samples the statistical trend is practically the same, which 
means that uncertainties in the velocity measurements due to inclination 
are not a major issue here. The black dots represent the total running 
median of our data and the error bars are determined by applying the 
bootstrapping statistical method. The different lines show the estimates 
in the literature for the total halo-to-stellar mass ratio once these have 
been scaled down by the factor specified in the labels. Lower frame: 
Mh/M,(< Rapi) as a function of Hubble type. The running median is 
overplotted in blue (error bars obtained via bootstrapping). 


thereby to the torque parameter Q\,. To do this accurately would 
require a sophisticated kinematic decomposition of observed ro¬ 
tation curves, but clearly such a study is not feasible given the 
size of the sample and the lack of rotation curve measurements 
for all of the galaxies. 


3.6.1. Halo-to-stellar mass ratio 


We obtained a first-order estimate of the halo-to-stellar mass 
ratio (Mh/Mt) inside the optical disk by comparing the circu- 
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lar velocity curve calculated from the 3.6 pm images with the 
inclination-corrected H i velocity amplitude (V^™). We assumed 
that at the optical radius Ropt (the radius enclosing 83% of the 
light in the blue band; 7?opt ~ 3.2 /ir for an exponential disk) the 
circular rotation velocity is close to the observed maximum ve¬ 
locity (y”“): 

(VSrf - Vl^^iRop,) + vLo(^opt). (12) 

From this we can infer 

^i6,™(^opl) 



MhlM,(< Ropt) 


■i 


where the factor F — 1.34 is the ratio between the mass con¬ 
tained by a spherical mass distribution and that enclosed by an 
exponential disk yielding a similar radial force at Ropt (Binney & 


Tremaine|1987|l. The gas contribution to the total rotation at the 


optical radius is here assumed to be negligible (e.g.|Rhee & van| 
Albadal 19961 |Verheijen| 19^) 1. 


In Fig. |6 we show the estimated halo-to-stellar mass ratio 
within Ropt as a function of stellar mass and compare this to 
various estimates in the literature for the total halo-to-stellar ra¬ 
tio. These estimates, based on abundance matching, weak lens- 
ing analysis in galaxy clusters, and halo occupation distribu¬ 
tion methods ([Moster et al.||2010t [Behroozi et aT]|2010| [Guo[ 
[et al.[[20T0 [Leauthaud et al.[ 2012[) , all predict a minimum of 
low-redshift Mh/M» at M» ~ IO'^^Mq. Our estimate for the 
halo mass within Ropt agrees well with the predictions if these 
are scaled down by a constant factor ~ 0.04. For the interval 
10'° Mq < M, < lO" Mq we obtain a mean Mh/M,(< Ropt) ~ 2 
and disks that are clearly less dominated by dark matter than the 
fainter counterparts (Mh/Mt(< Ropt) ~ 8 for M, < 10® Mq). 

The distribution of MhlMt{< Ropt) as a function of the re¬ 
vised Hubble stage is shown in the lower panel of Fig. Ex¬ 
treme late-type systems (Scd-Im) show a high dark matter rel¬ 
ative content within the optical disk (median of ~ 3.6 + 0.2), 
which is twice the median Mh/M»(< Ropt) of earlier-type sys¬ 
tems (~ 1.8 + 0.2 if 7° < 5). Although measurement uncertainties 
are larger among early-type SOs because of their gas-poor nature, 
lenticulars and early-type systems show a ~ 15% higher rela¬ 
tive dark matter content compared to intermediate-type galax¬ 
ies. This higher dark matter fraction among latest-type systems 


seems consistent with the reported trend in Falcon-Barroso et al. 
([2015[ see their Fig. 3) that was based on the CALIFA survey. 


3.6.2. Halo correction on the force profiles 

Figurej^indicates that the dark matter contribution inside the op¬ 
tical radius is not negligible, in particular for galaxies with stel¬ 
lar masses below IO'^Mq. Assuming that the halo distribution is 
spherically symmetric, the dilution of the normalized tangential 
force is 

gha1o-con(^) ^ . (14) 

RrCO + RhaloCO 


To study the influence of the dark halo in our bar strength mea¬ 
surements, we need to estimate its force profile inside the visible 
disk. In what follows we use several approximations, based on 
the so- called universal rotation curve models (URC) [Persic et al[] 


(jl996 PSS hereafter), combined with our S^G measurements 


and Hi amplitudes from literature. In Appendix [P] we compare 
in more detail how the URC fits the kinematic measurements for 
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NGC 4548 




Fig. 7. Ad hoc RC model for the barred galaxy NGC 4548 and the effect of the dark matter halo on the force profile. Left panel: Solid curves 
indicate the RC model based on the stellar component of the rotation curve, inferred from the 3.6 pm photometry. The amplitude of the URC 
halo component (isothermal sphere) has been corrected so that the total rotation curve matches the H i maximum velocity at the optical radius 
(Rop, = 3.2 /zr). The dotted line displays the URC model (exponential disk + isothermal sphere), while the dashed lines correspond to the models 
used in B2004. Central panel: Qj profiles, with and without halo correction. Right panel: Halo-corrected Qi profiles resulting from a measurement 
error in input parameters. A 20 % overestimate on the input parameters is assumed (in the mass-to-light ratio, the distance to the galaxy, and the 
observed velocity amplitude). We also assess the impact of reducing by one mag the /-band total luminosity (affecting the core radius in the URC 
models). The effect of matching the observed maximum velocity at 2.2/iR instead of at 3.2/2 r is tested as well. 


the sample galaxies, and we provide further details of our 
halo density profile determination. 

The URC models are based on a fit of exponential disk -H 
isothermal halo model on measured velocity profiles (see Ap- 
pendix[^for details). The isothermal halo implies 


(x^ + a^y ^ ^ 

where x - rlR^px is the radius normalized to the optical radius 
and Vca (halo velocity amplitude) and a (halo core radius) are 
given in terms of total luminosity in URC. Several corrections 
of the bar force accounting for the halo contribution are contem¬ 
plated; 

i) The first estimate for halo dilution on bar forces would be 
to replace the total radial force in the denominator of Eq. p^with 


the force predicted by the URC (Eq. D.l[ ). According to Eig. |D.l[ 
this should be more or less acceptable in a statistical sense. How¬ 
ever, it would completely discard the information about the indi¬ 
vidual stellar contribution to the velocity curve calculated from 
the S^G image. 

ii) A better estimate is to use the disk velocity curve as cal¬ 
culated from the 3.6 pm images and use URC solely to estimate 
the halo contribution. This was also the approach in |Buta et al.| 
( |2004| B2004 hereafter), where the halo core radius was taken 
from the URC model, based on the galaxy’s total blue band lu¬ 
minosity, while the halo amplitude Voa was obtained by fitting 
the relative contributions of the disk and halo at to those in 
the URC-models. Halo parameters were inferred from the URC 
halo-to-stellar mass relation. 

iii) However, compared to B2004, where no velocity mea¬ 
surements were available for the galaxies, we now have H i ve¬ 
locity amplitudes for most of the sample galaxies. We may use 
this by fitting Voa from the requirement that the modelled max¬ 
imum velocity within //opt matches the Hi velocity amplitude. 
In addition, we estimate the halo core from the /-band lumi¬ 


nosity based on the URC parameterization from Hendry et al. 
( 1997| l. Altogether, such an improved estimate can be made for 
484 galaxies with usable H i data. 


Our final halo models imply that ~ 10% of the disks in our 
sample are maximal according to the criterion of |Carignan &| 
|Ereeman| ( |1985| l and [van Albada et al.| ( |19851 l, based on the ve¬ 
locity of the stellar component relative to the total rotation eval¬ 
uated at 2.2/ir (for maximal disks Vdisk/^totai(2.2/jR) > 0.85). 

In the left panel of Eig. we show the RC model for the 
barred galaxy NGC 4548. This is based on the combination of 
stellar disk rotation curve combined with the estimated halo pro¬ 
file based on approximation (iii) above (see also Eig. |B.1| for 
an example with a more prominent dark matter halo). The mid¬ 
dle panel then indicates the corresponding halo-corrected force 
profile. Eor this particular early-type spiral galaxy, the effect of 
the halo is weak at rgb and the deviation from the non-corrected 
value is smaller than the uncertainty related to the unknown ver¬ 
tical thickness. However, for late-type galaxies, which typically 
host a higher relative amount of dark matter, the halo correction 
will have somewhat larger contributions, as discussed below in 
Sect. 5.3.1 In addition, in the right panel of Ei g. [7| we check 
the robustness of the halo model (see Appendix |E)|for a more 
detailed explanation) by assessing the uncertainties on the used 
input parameters, namely the mass-to-light ratio, the distance to 
the galaxy, and the observed velocity. In addition, we check the 
effect of matching the observed maximum velocity at a closer 
radius (2.2/ir instead of 3.2/ir). None of these changes result in 
a strong deviation in the (•< io%). 


4. Estimation of bar iengths and eiiipticities 

We used bar length measurements based on ellipse fitting (r^) 
and on visual inspection of the 3.6;um images (rvis), taken from 
HE2015. We also used and rqb values obtained from our 
force calculation. 

The analysis of ellipticity profiles in HE2015 was made for 
all those galaxies for which bars were identified in the morpho¬ 
logical classification of B2015 for the complete S^G sample. 
Based on P4 ellipse profiles, using the ellipticity maxima in the 
bar region and the constancy of the position angle in that region 
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NGC 4548 



-100 -50 0 50 100 


P4 Ellipticity radial profile 



P4 PA radial profile 



Radius (arcsec) 

Fig. 8. NGC 4548 bar characterization based on P4 isophotal ellipticity 
profiles. Top panel: Image of NGC 4548, in magnitude scale with range 
17-25 p 3 , (,fim(AB). The frame is aligned in a way that the y-axis points 
north and east is to the left. Axes are in units of arcsec. The ellipse 
corresponds to the maximum ellipticity isophote, while the straight line 
shows the visual estimate of the bar length and PA. Middle panel: ellip¬ 
ticity profile, with the maximum highlighted with a vertical line. Bottom 
panel: PA profile, with displayed with a vertical line. 


( jWozniak & Pierc^|1991[ |Wozniak et al.|[T995] l, 654 bars were 
identified in the 860 barred galaxies in B2015. Bar parameters in 
HE2015 were measured with two different methods: in the first 
method the bar length, position angle, and ellipticity were taken 
from the isophote fit corresponding to the ellipticity maximum in 
the bar region. In the second method, the length of the bar and its 
position angle were visually marked on top of the image. In both 
cases the bar was measured on the original sky image. The bar 
parameters (length, PA, e) were then converted to the disk plane, 
using the 2D analytical de-projection code in|Salo et al. \\999\ 
and LS2002 (the formulae are the same as in |Gadotti et al.|2007| l. 
Since the visual measurement of the bar only gave its length and 
position angle, the maximum of isophotal ellipticity was used in 
the conversion to the disk plane. 

An example for obtaining bar lengths from ellipse fitting is 
illustrated in Fig.|^for NGC 4548. The length is the radius of the 
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Fig. 9. Comparisons of different measurements of bar sizes, with the 
symbols colour-coded according to the different Hubble types. Top 
panel: Bar length calculated from the radius of maximum ellipticity and 
visual estimate (both converted to intrinsic disk plane values). Central 
panel: Visual bar length compared to the radius of maximum gravita¬ 
tional torque. Bottom panel: Visual bar length compared to the radius 
of maximum A 2 Fourier amplitude. 


maximum ellipticity, r^. As an upper limit of bar length, |Erwin| 
& Sparke ( |2003| ) and Erwin ( |2005[ E2005 hereafter) have used 


the first minimum in the ellipticity profile after r^, named as Lbai 
However, as discussed in |MicheTDansac & Wozniak| ( |2006| l, the 
location and value of such a minimum depends on the galaxy 
inclination and/or the real non-axisymmetric shape of the disk 
and the type of bar profiles, which is obvious also from the S"^G 
images. Therefore Lbai is not used in this study. If not otherwise 
mentioned, in the following the visually estimated bar lengths, 
after conversion to the disk plane, are used and referred as rbai - 
Furthermore, rbar always refers to the bar radius and not to the 
full length. 

Studies based on both observations (|Wozniak & Pierce|199ri 


Laurikainen & Salo 2002 Erwin & Sparke 2003|l and on N-bo^ 


simulations (Rautiainen & Salo|1999 


[Hs]|2002[ |Michel-Dansac & Wozniak 


jAthanassoula & Misirio-| 

|2006| l show that tends 


to underestimate the real size of bars. Using synthetic images, 
jAguerri et al.| ( |2009] l demonstrated that the length actually de¬ 
pends on the model describing the surface brightness profile of 
the bar (e.g. larger underestimation for Ferrers bars compared to 
Freeman and flat bars). 
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A statistical comparison between the bar sizes calculated vi¬ 
sually and the different independent systematic measurements 
of bar lengths (r^, rp^, and rQb) is presented in Fig. The 
measurements correlate tightly with We also confirm the 
tight correlation between the visual estimation of the bar size 


and both rpa ((Athanassoula & Misiriotis 2002 1 and rgb (Lau- 
[rikainen et al.||2i)D2| l. However, the later relation is much more 
scattered among the late-type galaxies. We obtain mean values 

< ^fi/^vis >= 0.97 (cr = 0.18), < rQb/rvis >= 0.62 (cr = 0.25) and 

< ^Az/Gis >= 0.75 (cr = 0.25). Thus, subject to natural measure¬ 
ment uncertainties, all the bar size proxies underestimate the bar 
length compared to the visual estimate, in agreement with what 
is reported in the literature. However, the difference between 
and Tvis is not as large as in previous studies. 

Finally, a small fraction of our galaxies show rpo. and rgb 
distances that fall beyond the bar length. Cases like this are 
seen for all Hubble types. These are cases in which the spiral 
arms strongly influence the potential in the bar region, where 
both maxima are determined by the combined contribution of the 
two non-axisymmetric stellar structures and occasionally reach 
beyond the bar limit. In the case of the ellipticity, bars with 
^e/Gis > 1 are explained by the presence of spiral arms and rings, 
as discussed in E2005 and|Comer6n et alI1(|2014|l. 


5. Characterization of bars as a function of Hubbie 
stage and famiiy 

In Table we present the mean values of the bar strengths (Q\,, 
e) and bar sizes (in physical units and normalized to R 25.5 
and /ir) for the different morphological classes and bar families. 

From now on, we refer to the galaxies in our sample as SOs if 
r < 0, early-type spirals if T e [0,3), intermediate-type spirals 
when T e [3,5), late-type spirals for T 6 [5,7], and Magellanics 
and irregulars for those with T > 1. With regard to the study of 
the bar properties in the Hubble sequence, we exclude ellipticals 
{T < -3) and also those galaxies that have a double morpholog¬ 
ical classification in B2015; this is the case of some lenticular 
and early-type spiral galaxies with an embedded disk-like inner 
structure with spiral(s) and a bar. However, there are only six 
such cases in our sample, and the measurements of bars are still 
provided, but are not included in any analysis. 

5. 1. Bar fraction in the Hubbie sequence 

We studied the bar fraction (/bar) for our sample of face-on and 
moderately inclined galaxies (i < 65°). Based on the meth¬ 
ods explained above, three different criteria were used to de¬ 
termine whether or not a galaxy has a bar: (1) from the visual 
classification in B2015, (2) by inspection of the radial ellipticity 
profiles, and (3) from the A 2 Fourier amplitude profiles and the 
presence of a well-defined four-quadrant butterfly pattern in the 
torque map. 

The upper panel of Fig. shows the histogram of all 
galaxies as a function of Hubble type and that of the barred 
galaxies, based on the three different criteria. It indicates that 
the sample is dominated by late-type systems, among which the 
frequency of bars is clearly decreased when criteria 2 and 3 are 
used instead of visual classifications. The lower panel shows the 
same, normalized to the total number of galaxies in the bin. 

A double-humped distribution of the bar fraction stands out 
regardless of the bar detection criterion, with a local minimum 
at r = 4. The physical reason for a lower bar fraction of these 
transitional Sb/c systems is not obvious, since they seem to be as 


All galaxies (2352) 
Inclination < 65°: 



-4 -2 0 2 4 6 8 10 



Fig. 10. Upper panel: Distribution of galaxies in our sample vs. their 
Hubble type. The black and orange dotted lines correspond to the his¬ 
togram of all the S'^G galaxies and those with inclination lower than 65°, 
respectively. For the galaxies with i < 65°, we display the histogram 
of barred galaxies in the Hubble sequence (solid lines), according to 
the various criteria explained in the text. Lower panel: Bar fraction vs. 
Hubble stage for a sample of face-on and moderately inclined galaxies 
(i < 65°). The fraction of galaxies with an identifiable maximum in the 
amplitude of their higher order even Fourier components is also shown. 
Uncertainties are calculated assuming a binomial distribution. 


massive and have similar dark matter fractions as earlier types of 
spiral galaxies in the S^G sample. Interestingly, no similar dip is 
present when the bar fraction as a function of M* is studied (see 
Sect.[54Tl. 

For later types, the bar fraction rises again. However, it ap¬ 
pears that more bars are identified visually than with the other 
methods. Indeed, the bar fraction is reduced by ~ 30% and 
~ 50% for r > 4 if e or A 2 and Qj profiles are used (criteria 
2 and 3 above, respectively). This is natural, taking into account 
that very weak bars that are still visible in the images can be 
lost in ellipse fitting analysis where clumpy star-forming regions 
can induce stronger ellipticity maxima than the bars. Similarly, 
dust lanes can distort the shape of the isophotes. Dust lanes and 
H II regions can also make A2 profiles noisy. Many of these late- 
type bars would possibly be overlooked if they were observed at 
higher redshifts, given their faint disks. 

The bar fraction is remarkably high for Sab (T - 2) galaxies, 
and there is a drop in the bar fraction at Hubble types earlier than 
T = 0, consistent with|Laurikainen et al.|(|20l3|). 
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Bar fraction: 

Visual 

Ellipse fitting 

3A““&eb 

Motal 

T 6 [-3,0) 

49.1 ±4.6% 

40.5 ± 4.6 % 

43.1 ±4.6% 

116 

T e [0,3) 

69.8 ± 3.2 % 

61.8 ±3.3% 

61.8 ±3.3% 

212 

T e [3,5) 

55.8 ± 3.6 % 

48.9 ± 3.6 % 

46.8 ± 3.6% 

190 

T e [5,7] 

75.1 ±2.4% 

60.2 ± 2.7 % 

50.8 ± 2.8% 

329 

T>1 

60.6 ± 2.2 % 

37.1 ± 2.2 % 

25.5 ± 2.0% 

498 


Table 1. The mean bar fractions based on different bar identification 
criteria and for different morphological types. Binomial counting errors 
and the total number of galaxies within the bins are indicated. 


Ohta et aLl ( 1990| l showed that for early-type galaxies, m - 4 
and m = 6 components are not negligible in the azimuthal den¬ 
sity profiles of barred galaxies. In the lower panel of Fig. 10 we 


also show the bars for which the maximum of the higher order 
component amplitudes could be detected in the bar region. Of 
the early-type systems (T < 4) for which and were cal¬ 
culated, approximately one half of the galaxies present reliably 
identifiable A™“, A™“ and A™“. In contrast, very few galaxies 
with r > 4 have a significant contribution from the higher-order 
even components to the Fourier decompositions of bars. 

The bar fractions based on different bar detection criteria and 
for different morphological types are listed in Table For com¬ 
plementary details about the visual bar fraction in the S'^G, see 
B2015. 

5.2. Bar lengths in the Hubble sequence 

Bar lengths are displayed as a function of the galaxy Hubble 
stage in Fig.[TT] In the upper panel the distribution of bar lengths 
in kpc is displayed, and in the second and third panels we nor¬ 
malized them by the galaxy size (as measured by R 25.5 and /zr). 
We confirm the result from |Elmegreen & Elmegreen| ( |1985] l, who 
showed that the bars in early-type spiral galaxies (T = 0-2) 
are typically longer than the bars hosted by late-type spirals 
(r = 3 - 7). 

More specifically, we observe that bar sizes tend to increase 
from r = 5 towards the Sa and SO/a galaxies (in agreement with 
e.g. Martini 1995 1 , and then drop among the SOs (as found ear¬ 
lier in E2005 and |Laurikainen et al. [20071 1. This is roughly the 
same for the absolute and relative bar sizes and is independent of 
the normalization used (either by R 25.5 or /zr), but the trend flat¬ 
tens among SOs when normalizing by the disk scale length. The 
difference between the disk-relative bar sizes of early-type and 
late-type spirals is not as pronounced as the contrasts of a factor 
~ 3 and ~ 2.5 reported in |Martin| ( |1995| l and E2005, respec¬ 
tively. Indeed, the strongest difference is found by comparing Sa 
to Sc galaxies, the former hosting bars that are ~ 1.75 larger 
on average. On the other hand, Sd galaxies also typically have 
longer bars than Sc galaxies. The mean values for the different 
morphological types and bar classes are summarized in Table 

The distribution of bar lengths in the Hubble sequence that 
we obtain is more or less consistent with the statistics in E2005 
(we have a sample larger by a factor ~ 6 ), who found a mean 
bar size of ~ 3.3 kpc for SO-Sb galaxies (with the bar sizes dis¬ 
tributed in a similar range as we observe, although we find bars 
shorter than 1 kpc among early-type galaxies). Nevertheless, his 
late-type galaxies (Sc-Sd) had bars with a mean size of ~ 1.5 
kpc (ranging from 0.5 to 3.5 kpc), while in their T range ours 
are typically ~ 2.5 kpc, spanning in a wider interval (0.5 - 6 
kpc). Altogether, our estimates give a mean bar size of ~ 2.5 kpc, 
early-type systems (S0“-Sb) hosting bars that are a factor ~ 1.5 
longer in physical size than those in late-type galaxies (Sbc-Im). 
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Fig. 11 . Bar length distribution in terms of the galaxy Hubble stage. 
The blue line corresponds to the running mean, with the standard error 
of the mean indicated with an error bar. The green vertical line delim¬ 
its the region of SO galaxies. The first row displays the bar length (in 
units of kpc) vs. T. The red line shows the mnning mean of the radius 
of maximum ellipticity, with an offset of +0.25 in the x-axis. Shown in 
the second and third rows are the distributions of bar length (visual es¬ 
timate), normalized by f?25.5 and /zr. To avoid overlapping of points, the 
T values in the j:-axis have been randomly displaced from the integer. 


Eor scaled bar sizes, the measurements in E2005 are made 
in intervals 0.2 - 0. 8/^25 and 0.5 - 2.5/zr for SO-Sab galaxies, 
with mean values ~ 0.38/^25 and 1.4/zr. In their T range we also 
observe shorter bars (as short as ~ 0 . 1 /? 25.5 and a few cases with 
^bai /ZiR < 0.5), and we obtain lower mean bar sizes (~ 0.27/?25.5 
and ~ 1.25/zr). The disk-relative bar sizes of late-type systems 
in E2005 are found in the ranges 0.05 - 0.35f?25 and 0.2 - 1 .5/zr, 
with mean values of ~ 0.14/?25 and ~ 0.6/zr. We also measure 
bars among Sc-Sd galaxies that can be as large as ~ 0.67? 25.5 and 
~ 2/zr, and in general we observe slightly longer bars in these 
systems than E2005. 

Aguerri et al. (|20051l found a mean relative bar length 


(rbai//zR) = 1.21 + 0.08 (as measured from the photometric 
decomposition including lenses). In addition, numerical simu¬ 
lations made by O’Neill & PubinsE] ( |2003[ ) and Valenzuela & 
|Klypin| ( |2003| l predicted (rbar/fiR) = 11 ± 0.33 and 1.0 - 1.2, 
respectively. In this work, SOs and early-type spirals are the only 
systems with < rbar/^R >> 1, while for galaxies with T > 3 
bar sizes relative to /zr are on average below the predicted val- 
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r <0 

T s [0,3) 

T e [3,5) 

T e [5,7] 

T s (7,10] 

^Qb/^25.5 

fQblhR 

0.18±0.01 

0.87±0.06 

0.21±0.01 

0.95±0.03 

0.17±0.01 

0.77±0.04 

0.09±0.00 

0.35±0.02 

0.13±0.01 

0.43±0.02 

fA21^25.5 
rA2lhR 

0.20±0.01 

0.97±0.06 

0.22±0.01 

0.99±0.04 

0.17±0.01 

0.79±0.04 

0.12±0.00 

0.47±0.02 

0.19±0.01 

0.70±0.04 

f €1^25.5 
rJhR 

0.22±0.01 

1.07±0.06 

0.26±0.01 

1.21±0.04 

0.21±0.01 

1.00±0.04 

0.19±0.01 

0.72±0.02 

0.27±0.01 

0.95±0.03 

^Qb vis 

0.75±0.03 

0.75±0.02 

0.77±0.02 

0.49±0.02 

0.50±0.02 

rtalr^h 

0.83±0.02 

0.79±0.02 

0.80±0.02 

0.66±0.02 

0.74±0.03 

f e/^vis 

0.91±0.01 

0.98±0.01 

1.03±0.02 

0.96±0.02 

0.95±0.01 


Table 2. Statistics of the scaled bar sizes as measured from rgb, rA 2 
and Tf, and comparison of the visual measurements with the rest of the 
estimates. It is presented the mean and standard deviation of the mean. 


ues. This is natural considering that early-type bars resemble the 
type of bars resulting from numerical simulations, as stated in 
E2005. 

The lower mass systems at the end of the Hubble sequence 
are interesting galaxies. As expected, the absolute bar sizes are 
smaller than in the bright galaxies (SO -She galaxies have on 
average the same size), but the normalized sizes are on average 
even larger (compared to intermediate-type spirals). In fact, there 
is a tendency of increasing bar length from T - 6 towards the 
irregular galaxies. 

It seems that among the late-type systems, the more clumpy 
the galaxy, the longer the relative size of the elongated struc¬ 
ture embedded in the underlying disk. However, it is worth to 
mention that even among the smallest galaxies, the bars are long 
enough (97% of them are longer than 10 pixels, with the shortest 
bar being 6 pixels long) so that they are not likely to arise solely 
from the joining of a pair or group of star-forming clumps, that 
is, they are not visual artifacts (as can happen with measurements 
of bars at high redshifts). This trend of increased relative size 
among the faint galaxies, which was also present in the observa¬ 
tions of Laurikainen et al.| ( |2007| l in spite of the poorer sampling 
of T-types > 7, is maintained regardless of the bar size estimate 
that we use (see also Table [^. In addition, we have divided our 
sample into two bins in distance and checked that the trends for 
the bar sizes are exactly the same, confirming that the angular 
resolution is not a problem here. 


5.3. e, Q^, andA 2 in the Hubble sequence 

The prominence of a bar is displayed as a function of Hubble 


type in Fig. [T2| In addition to Qb and discussed above, we 


also used the ellipticity maximum in the bar region, e, as a proxy 
of the bar strength. 

We conhrm the tendencies from the previous studies, which 
showed that for the bright spiral galaxies with Hubble types 
T=0-6, the mean Qb increases towards the later types, while 
the mean decre ases (likewise A™“). Our results agree with 
those of |Buta et al.| ( |20T0| l, who showed that the frequency of 
strong bars, as estimated with Qb, is lower in SO systems. Early- 
type SOs are also characterized by a lower A™“. On the other 
hand, the maximum ellipticity is maintained nearly constant in 
each Hubble type bin, in agreement with Marinova & Jogee 
( |2007| l. Here we additionally show that e decreases from T=0 
towards the early-type SOs, as found earlier inlLaurikainen et al.l 
( |2007l l. 

The behaviour of the torque parameter among early-type 
galaxies is mostly explained by the dilution of Qb by the under¬ 
lying radial force field, primarily generated by bulges (so-called 
bulge dilution. Block et al.|2001[ LS2002). The high Qb values 
in the low-mass galaxies can be understood for the weak un¬ 


derlying disks (hence low Er, which is the denominator in the 
equation of Qb). Furthermore, we have already given additional 
evidence that dark halos dominate the mass distribution of these 
systems, whose effect on Qb is discussed in Sect. |5.3.T| 

We further investigate the loci of rA 2 , rgb , and with respect 
to the visual bar length in the lower three panels of Fig. 12 and 
also in Tablel^with galaxies separated in morphological classes 
(see also Fig.M. As we have already discussed, these three prox¬ 
ies of the bar size tend to underestimate the bar length as com¬ 
pared to the visual measurement. In addition, both rAilr^is and 
rgb/Gis ratios drop for late-type galaxies. This drop is remark¬ 
ably large for the maximum gravitational torque and is related 
to the reduced central concentration towards later-type systems, 
causing rQb to move inwards (see ILaurikainen et al.||2002|). In 
spite of this, all of these bar length proxies have the same be¬ 
haviour in the Hubble sequence. 

How much the above tendencies depend on the mass of the 
galaxy or on the masses of the different structure components in 
galaxies is discussed in Sect.|^ 


5.3.1. Sources of uncertainty on the Qb estimate 

We have carried out tests to probe possible sources of uncer¬ 
tainty in our Qb estimates; (1) the assumption of a constant scale 
height, (2) the impact of non-stellar emission in the bar region, 
(3) the contribution of the spiral arms torques to the overall Qb 
value, and (4) the effect of dark matter halos on the bar forcing. 

Here, we check whether these uncertainties can be large 
enough to change the observed trend of the bar torque param¬ 
eter in the Hubble sequence (see Fig.[T3|). 

We have already concluded that uncertainties associated to 
the constant assumption and the effect of non-stellar contam¬ 
inants are similar in size as the uncertainty related to the thick¬ 
ness of the disk (~ 10 - 15%). However, as the effect of the 
dark matter halo on the bar force depends on whether our disks 
are maximal or sub-maximal, the halo correction is likely to in¬ 
crease with T. This is indeed the case for T > 4. Interestingly, 
earlier-type systems also have a strong halo correction, which 
is fairly constant in the T-type range [-3,4]; the amplitude of 
the halo correction on Qj increases with radius and the loci of 
rQb with respect to the galaxy centre moves outwards for more 
concentrated galaxies, explaining the lack of difference on the 
mean gh^io-corr between intermediate-type spirals and early-type 

spirals and SOs. We obtain (e“°^“"'/eb) = 0.83 (cr = 0.15). 

Here, we also compare the approaches described in 
Sect. |3.6.^ for the halo correction. As shown in Fig.[T^ the halo 
correction made in B2004 (URC halo + circular velocity from 
the potential) and the direct correction using the URC halo-i-disk 
yield quite similar mean values. Compared to the improved 
method applied in this work, we obtain a slightly larger correc¬ 
tion for the earliest types. The smaller Qb correction present in 
the data of B2004 (6%) might depend on the use of a sample 
with brighter galaxies and also on the methodology. 

In the statistical sense the bar-only force measurements do 
not deviate considerably from the raw Qb measurements (~ 10% 
on average), in agreement with |Laurikainen et al.| ( |2007] l. Only 
Sb-Sc galaxies show values that are slightly higher than 

the uncertainty in the raw estimate of Qb. Setting the m > 0 
Fourier components to zero beyond 1.2 times the bar size for the 
potential calculation (instead of using rem = rb^Q gives practi¬ 
cally the same results as using the raw Qb. 

In conclusion, none of the analysed sources of uncertainty 
alter the statistical trend of Qb in the Hubble sequence. The dif- 
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Fig. 12. Different bar strength indicators as a function of the Hubble Type. Upper row: The 2b, A““, A™’' and e are plotted vs. the integer value 
of the revised numerical Hubble stage (from left to right, respectively). Lower row: Radius of maximum gravitational torque, Fourier amplitude 
and ellipticity, normalized to the visual estimate of the bar length as a function of the Hubble stage. The running mean is overplotted with a solid 
line and the standard error of the mean is indicated with a vertical error bar. T values in the x-axis (integers) have been randomly displaced for the 
sake of avoiding point overlapping. For all the plots, the green vertical line demarcates the region of SO galaxies. In the central upper panel, the 
moving average of bar maximum amplitude of the m = 4 Fourier component (A““) is also displayed. 


The vertical thickness of the disk 



-4 -2 0 2 4 6 8 10 

T 


Dark matter correction Bulges and Spiral Arms 



Fig. 13. Running mean of 2b in the Hubble sequence, with and without spiral arms and dark halo corrections. Left panel: Mean bar torque 
parameter obtained using the two disk thickness estimates applied in this study (see Sect. |3.1.2^ : [de Grijs|(T998) and |Speltincx et ar] ( |2008| l (blue 
and red solid lines, respectively.). The lower and upper limits for the disk thickness are calculated using the de Grijs’ relation, from which the 
corresponding mean bar torque parameter delimits the uncertainty in the estimate of 2b (blue shaded area). Central panel: The mean torque after 
the halo correction, following the different approaches described in Sect. |3.6.^ (i) using the universal rotation curve, (ii) following the methodology 
of B2004, and (iii) adding V '3 b^n]+URC halo, and correcting the URC halo amplitude to match the observed velocities. The unrealistic extreme 
case (iv) of assuming Fbaio/F^6^,^ to be constant at all radii is also shown. Corrections (i), (ii), (iii), and (iv) correspond to the light green, orange, 
red, and dark green solid lines, respectively. The shaded area is computed from all the reliable bar strength measurements presented in this paper, 
while the remaining moving averages (including that one of the raw 2b) are computed using only the galaxies with available halo correction. Right 
panel: The moving average of 2b after the spiral arm correction. Galaxies with prominent bulges are studied separately, showing the role they play 
in diluting the bar forces. We obtain setting the m> 0 Fourier components to zero (1) beyond the bar length and (2) beyond 1.2 times the 

bar length in the potential calculation. 
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Fig. 14. Moving average of the gravitational force evaluated at the end 
of the bar - QxCfbar) - with and without halo correction. 



SB SAB SAB SAB 


Bar classes 

Fig. 16. Mean bar strengths and disk-relative bar sizes for all the differ¬ 
ent bar families. Error bars are the standard deviation of the mean. Both 
irregular and spiral galaxies are considered in the analysis. 


ference of a factor ~ 2 in Qb between early- and late-type sys¬ 
tems is slightly reduced when taking into account the dark halo 
correction, but this is not sufficient to change the monotonically 
raising trend of Qb in the Hubble sequence. In spite of this, for 
individual galaxies dark halos might change the bar force mea¬ 
surements substantially. 


5.3.2. Gravitational torque at the end of the bar 


In Fig. we assess how strong the halo correction would be 
if we were to evaluate the gravitational force at a radius larger 
than rgb. For instance, this is of vital importance to shed light 
on the bar-spiral interplay based on comparisons between their 
strengths, an idea which has remained controversial (examples 
either supporting or refuting a connection in the coupling of 
these stractures can be found in e.g. Seigar et al.||2003t |Block| 


|et al.|20()4{ |Buta et al.|2005[ |Durbala et al.|20d^ although the 
recent analysis carried out in Buta et al. (20091 and jSalo et aT7 


( | 2010 | l found evidence of a correlation between local bar forcing 
and local spiral amplitude. 

We observe that the behaviour of SrC^bai) in the Hubble se¬ 


quence is similar to the distribution of Qb shown in Fig. 14 with 
certain flattening among the later-type systems. This is explained 
by the aforementioned tendency of rgb to move inwards for the 
faintest disks, which causes Qj to experience a more pronounced 
drop at rbai relative to 2 b- 


The halo correction seems to be of the same order as in Qb 
for T < 5 and larger for later types (~ 30 - 50%), making 
the average 2 !j?‘'°“‘^™(rbar) value roughly constant (~ 0 . 2 ) for 
intermediate-type spirals and later-type systems, and following 
a trend that is fairly similar to the trend of the bar intrinsic ellip- 
ticity. 


5.4. Bar parameters and the family class 

While comparing the visual bar classihcations in B2015 (using 
the de Vauculeurs family classihcation nomenclature from RC3) 
with the measures of the prominence of the bar obtained in this 
study, it is evident that the quantitative estimates are well con¬ 
sistent with the visual ones. This is shown in Figs. [TS] and [Tb] 


(see also Table |^, where the distributions of bar strengths and 
lengths are shown separated for the different bar families. 


Erwin I (|2005[) concluded that strong (SB) and weak (SAB) 


bars hosted by early-type galaxies differ primarily in ellipticity, 
being very similar in size, while SB bars in late-type spiral galax¬ 
ies are twice as large as SAB bars. We observe that SB, SAB, 
SAB, and SAB bars differ in ellipticity as well as in and 
for all morphological types, as shown in Table The decreas¬ 
ing tendency from SB towards the SAB family is particularly 
clear for the bar torque parameter. In the statistical sense, there is 
also a difference between weak and strong bars for the bar sizes, 
as visually categorized; SB and SAB galaxies are ~ 10 - 15% 
longer relative to the disk size. We already hnd a similar dif¬ 
ference in disk-relative bar sizes between weak and strong bars 
among early-type spirals (0 < T < 3). However, for SOs weak 
and strong bars have similar mean bar sizes. On the other hand, 
although our bars are on average slightly longer for SB and SAB 
families in the late-type systems, we do not obtain such a huge 
difference between weak and strong bars (factor ~ 2 in E2005). 


5.5. Comparison of bar strength measurements 

In Eig. [T^ we show a comparison of the different bar strength 
indexes, conhrming some of the trends earlier reported in the 
literature. We observe a correlation between the intrinsic ellip- 
ticity of the bar and the bar gravitational torque (e.g. [Block et al. 
[2001 1 LS2002), which shows that underlying potential tends to 
adjust the orbit of the stars that constitute the bar. Galaxies of all 
Hubble types present the same pattern, although two families of 
bars can be identihed based on the trends in the plot: SOs, early- 
type, and intermediate-type spirals occupy the lower part, and 
late-type spirals, Magellanics, and irregulars are located in the 
upper area. Eor these families the data are htted separately with 
a polynomial of order 2 : 

2b = 0.1 - 0.34e +l.\6s^ , if T < 5 and s 6 (0.15,0.8), 

2b =0.21-0.49e-H 1.46 s2 , if T > 5 and e 6 (0.25,0.9). 


Eor the galaxies with T <5, the deviation with respect to the 
htting curve is small (~ 5%), while for later types the scatter is 
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Table 3. Statistics of the measured bar strengths and lengths (in kpc and normalized to /?25 5 and / jr ), for all barred galaxies and also separated into 
morphological types and families. For each of the different measurements and subsamples, the mean value, the standard deviation of the mean and 
the number of galaxies within the bin are given. 


Buta et al. 2015 

Qb 

p^max 

€ 

rbar/R25.5 

Ibar/hR 

rbar(kpc) 

Barred 

0.35 ± 0.01 (576) 

0.41 ± 0.01 (587) 

0.52 ± 0.01 (654) 

0.24 ± 0.00 (7S2) 

0.98 ± 0.01 (742) 

2.55 ± 0.06 (784) 




Morphological type 




r <0 

0.16 ±0.01 (50) 

0.45 ± 0.03 (53) 

0.40 ± 0.02 (47) 

0.24 ± 0.01 (56) 

1.22 ±0.06 (50) 

2.91 ± 0.25 (56) 

T 6 [0,3) 

0.26 ± 0.01 (132) 

0.51 ±0.02 (134) 

0.51 ±0.01 (131) 

0.27 ± 0.01 (746) 

1.24 ±0.03 (137) 

3.76 ±0.19 (746) 

T s [3,5) 

0.28 ± 0.01 (89) 

0.39 ± 0.02 (95) 

0.52 ± 0.02 (93) 

0.20 ± 0.01 (105) 

0.96 ± 0.03 (702) 

2.87 ±0.18 (705) 

T 6 [5,7] 

0.42 ± 0.02 (168) 

0.36 ± 0.01 (174) 

0.53 ± 0.01 (198) 

0.20 ± 0.01 (234) 

0.76 ± 0.02 (230) 

2.07 ± 0.07 (255) 

T >1 

0.45 ± 0.01 (137) 

0.37 ±0.01 (131) 

0.54 ± 0.01 (185) 

0.28 ± 0.01 (247) 

0.99 ± 0.02 (223) 

2.06 ± 0.07 (242) 

SO- 

0.10 ±0.02 (10) 

0.33 ±0.04 (77) 

0.24 ± 0.03 (7) 

0.21 ± 0.02 (74) 

1.13 ±0.12 (72) 

2.61 ± 0.58 (74) 

SO® 

0.16 ±0.03 (16) 

0.41 ± 0.04 (18) 

0.39 ± 0.03 (76) 

0.23 ± 0.02 (18) 

1.30 ±0.10 (77) 

2.28 ± 0.42 (18) 

S0+ 

0.17 ±0.01 (24) 

0.54 ± 0.05 (24) 

0.45 ± 0.02 (24) 

0.27 ± 0.01 (24) 

1.21 ±0.09 (27) 

3.55 ± 0.32 (24) 

so/a 

0.23 ± 0.02 (40) 

0.54 ±0.04 (47) 

0.48 ± 0.02 (42) 

0.26 ± 0.02 (47) 

1.21 ±0.06(45) 

3.88 ± 0.35 (47) 

Sa 

0.26 ± 0.02 (47) 

0.52 ± 0.04 (49) 

0.49 ± 0.02 (45) 

0.28 ± 0.02 (54) 

1.27 ±0.07 (49) 

3.78 ± 0.32 (54) 

Sab 

0.27 ± 0.02 (41) 

0.44 ± 0.03 (40) 

0.54 ± 0.02 (40) 

0.26 ±0.01 (47) 

1.24 ±0.05 (47) 

3.47 ± 0.30 (47) 

Sb 

0.30 ± 0.02 (50) 

0.44 ± 0.03 (54) 

0.55 ± 0.02 (52) 

0.22 ± 0.01 (58) 

1.05 ±0.05 (56) 

3.13 ±0.24 (58) 

She 

0.26 ± 0.02 (36) 

0.31 ± 0.03 (38) 

0.49 ± 0.03 (38) 

0.18 ±0.01 (43) 

0.84 ± 0.05 (42) 

2.50 ± 0.28 (45) 

Sc 

0.35 ±0.03 (51) 

0.32 ± 0.02 (52) 

0.49 ±0.02 (67) 

0.16 ±0.01 (68) 

0.69 ± 0.03 (67) 

1.84 ± 0.10 (6S) 

Sed 

0.43 ± 0.02 (59) 

0.37 ±0.02 (61) 

0.53 ± 0.02 (63) 

0.18 ±0.01 (74) 

0.72 ± 0.03 (73) 

2.06 ±0.12 (75) 

Sd 

0.48 ± 0.03 (57) 

0.39 ± 0.03 (60) 

0.55 ± 0.02 (73) 

0.24 ±0.01 (97) 

0.86 ± 0.04 (89) 

2.26 ±0.13 (97) 

Sdm 

0.47 ± 0.02 (74) 

0.39 ± 0.02 (70) 

0.57 ± 0.02 (95) 

0.27 ± 0.01 (777) 

0.93 ± 0.03 (114) 

2.30 ± 0.10 (77S) 

Sm 

0.43 ± 0.02 (49) 

0.38 ± 0.02 (47) 

0.52 ± 0.02 (64) 

0.29 ± 0.01 (84) 

1.02 ±0.05 (72) 

1.95 ± 0.12 (S4) 

Im 

0.38 ± 0.04 (12) 

0.28 ± 0.05 (72) 

0.48 ± 0.04 (24) 

0.30 ± 0.02 (35) 

1.11 ±0.06(54) 

1.57 ± 0.15 (55) 




Galaxy family 




SB 

0.45 ± 0.01 (260) 

0.48 ±0.01 (267) 

0.59 ± 0.01 (277) 

0.25 ± 0.01 (323) 

1.02 ±0.02 (308) 

2.73 ±0.10 (525) 

SAB 

0.34 ± 0.02 (67) 

0.39 ± 0.02 (69) 

0.55 ± 0.02 (76) 

0.27 ± 0.01 (96) 

1.03 ±0.04 (89) 

2.38 ± 1.44 (96) 

SAB 

0.26 ± 0.01 (193) 

0.38 ± 0.02 (193) 

0.46 ± 0.01 (226) 

0.24 ± 0.01 (264) 

0.97 ± 0.02 (255) 

2.54 ±0.10 (266) 

SAB 

0.20 ± 0.02 (33) 

0.26 ± 0.02 (37) 

0.36 ± 0.02 (46) 

0.18 ±0.01 (63) 

0.78 ±0.06 (67) 

2.06 ±0.19 (65) 




SO galaxies: T < 0 




SB 

0.19 ±0.02 (17) 

0.50 ± 0.04 (77) 

0.44 ± 0.02 (18) 

0.25 ± 0.02 (18) 

1.31 ±0.07 (77) 

2.88 ± 0.43 (18) 

SAB 

0.17 ±0.03 (7) 

0.42 ± 0.05 (8) 

0.39 ± 0.04 (7) 

0.18 ±0.02 (9) 

1.13 ±0.10 (8) 

2.14 ± 1.42 (9) 

SAB 

0.13 ±0.01 (17) 

0.46 ± 0.07 (77) 

0.39 ± 0.03 (77) 

0.24 ± 0.02 (18) 

1.10 ±0.11 (76) 

2.79 ± 0.39 (18) 

SAB 

0.13 ±0.03 (8) 

0.39 ± 0.07 (8) 

0.32 ± 0.09 (4) 

0.29 ± 0.03 (8) 

1.45 ±0.23 (7) 

3.98 ± 1.00 (S) 



Early-type spirals [SO/a-Sb): 

T e [0,3) 



SB 

0.33 ± 0.02 (54) 

0.59 ± 0.03 (55) 

0.59 ±0.01 (57) 

0.29 ± 0.01 (59) 

1.31 ±0.05 (56) 

4.15 ±0.31 (59) 

SAB 

0.28 ± 0.05 (8) 

0.57 ± 0.08 (8) 

0.56 ± 0.03 (8) 

0.30 ± 0.02 (8) 

1.36 ± 0.14 (S) 

4.23 ± 0.86 (8) 

SAB 

0.20 ± 0.01 (58) 

0.47 ± 0.03 (57) 

0.46 ± 0.02 (56) 

0.26 ±0.01 (67) 

1.22 ±0.05 (56) 

3.63 ± 0.30 (67) 

SAB 

0.19 ±0.08 (3) 

0.23 ± 0.07 (4) 

0.34 ± 0.09 (5) 

0.23 ± 0.03 (7) 

1.05 ± 0.17 (7) 

2.90 ± 0.62 (7) 



Intermediate-type spirals [Sb-Sc): T 6 [3,5) 



SB 

0.43 ± 0.03 (39) 

0.49 ± 0.03 (38) 

0.61 ± 0.02 (39) 

0.23 ±0.01 (41) 

1.02 ±0.05 (47) 

3.30 ± 0.36 (47) 

SAB 

0.28 ± 0.02 (17) 

0.37 ± 0.05 (77) 

0.58 ± 0.03 (76) 

0.22 ± 0.02 (17) 

1.02 ±0.08 (77) 

2.88 ± 1.96 (77) 

SAB 

0.27 ± 0.02 (59) 

0.36 ± 0.02 (64) 

0.49 ± 0.02 (67) 

0.18 ±0.01 (78) 

0.83 ± 0.03 (77) 

2.29 ± 0.13 (7S) 

SAB 

0.23 ± 0.04 (14) 

0.25 ± 0.02 (76) 

0.39 ± 0.04 (79) 

0.12 ±0.01 (23) 

0.63 ± 0.07 (25) 

1.59 ±0.12 (25) 



Late-type spirals [Sc-Sd]: T 

6 [5,7] 



SB 

0.52 ± 0.02 (92) 

0.44 ± 0.02 (93) 

0.61 ± 0.01 (94) 

0.21 ± 0.01 (108) 

0.80 ± 0.03 (107) 

2.26 ±0.12 (70S) 

SAB 

0.41 ±0.04 (12) 

0.34 ± 0.02 (13) 

0.62 ± 0.02 (75) 

0.22 ± 0.03 (75) 

0.72 ± 0.06 (74) 

2.14 ±0.91 (75) 

SAB 

0.32 ± 0.02 (46) 

0.29 ± 0.02 (49) 

0.47 ± 0.02 (60) 

0.20 ± 0.01 (74) 

0.79 ± 0.04 (75) 

1.97 ± 0.11 (75) 

SAB 

0.25 ± 0.06 (8) 

0.17 ±0.02 (9) 

0.32 ± 0.03 (18) 

0.13 ±0.01 (24) 

0.60 ± 0.06 (24) 

1.46 ±0.12 (24) 



Magellanic and 

Irregular galaxies [Sdm,Sm,Im]: T >1 



(S/I)B 

0.52 ± 0.02 (74) 

0.42 ± 0.02 (73) 

0.60 ± 0.02 (89) 

0.27 ± 0.01 (113) 

1.00 ± 0.04 (705) 

2.12 ±0.10 (775) 

(S/1)AB 

0.42 ± 0.02 (25) 

0.37 ± 0.03 (25) 

0.54 ± 0.03 (32) 

0.31 ±0.01 (49) 

1.04 ±0.05 (44) 

2.02 ± 1.16(49) 

(S/1)AB 

0.35 ± 0.01 (33) 

0.29 ± 0.02 (28) 

0.47 ± 0.02 (52) 

0.30 ± 0.01 (62) 

0.99 ± 0.04 (60) 

2.05 ±0.13 (65) 

(S/1)AB 

0.28 ± 0.04 (4) 

0.19 ±0.02 (4) 

0.39 ± 0.06 (9) 

0.23 ± 0.02 (72) 

0.70 ±0.08 (77) 

1.70 ± 0.23 (72) 
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Fig. 15. From left to right and top to bottom, histograms of the distributions of the bar gravitational torque parameter Q^,, bar density amplitude 
A™", bar intrinsic ellipticity s, f? 25 . 5 -relative bar size, /iR-relative bar size, and bar length in physical units. 
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Fig. 17. Comparison of bar strength indexes. Galaxies are separated into T bins. Spearman’s rank correlation coefficients (together with the 
significance) are indicated in all the plots for the total sample (see Tablej^for the different morphological types). Systems with relatively massive 
(faint) bulges are plotted with filled (empty) circles. In the left and central panels, the dotted and dashed lines display a second-order polynomial 
fit to the data separately for F <5 and T > 5, respectively. In the right panel we show the linear fits for different morphological classes. 


Correlations; 

Qb vs. £ 

£ VS. 

gb vs. A™ 

gnmo-corr ^ 

Ct(''.) vs. e 

vs. £ 

g.„-c„„vs.A,"" 

fiT(''bT) vs. A 2 (rb„) 

er"-“""(rb„) vs. A 2 (rb„) 

r <0 

0.80(6.5- 10-") 

0.84(5.2- 10-") 

0.72 (3.9 - 10-“) 

0.79(8.1 - 10-“) 

0.78(9.0- 10-") 

0.78 (9.0 - 10-") 

0.73 (3.2 - 10-’) 

0.40 (0. - 10-") 

0.44(0- 10-") 

T e [0,3) 

0.88 (3.1 ■ 10-*') 

0.58(1.3 - lO-"") 

0.62 (3.4 - 10-‘“) 

0.80(3.2- 10-"") 

0.85 (5.9 - 10-"') 

0.85 (5.9 - 10-"') 

0.64(9.8- 10-'") 

0.43 (3.7 - 10-*) 

0.42 (4.7 - 10-“) 

T e [3,5) 

0.88(1.6-10-^') 

0.65(1.9-10-") 

0.70 (3.9 - lO-") 

O 

oq 

00 

o 

0.68 (7.3 - 10-") 

0.68 (7.3 - 10-") 

0.67(2.4- 10-") 

0.42 (6.3 - 10-“) 

0.44(1.3- 10-’) 

T e [5,7] 

0.86 (7.5 - 10-") 

0.81(7.9-10-"“) 

0.88 (0.0) 

0,81 (3.0- 10-^^) 

0.70(1.2- 10-"“) 

0.70(1.2- 10-"“) 

0.83(1.0- 10-"") 

0.58 (8.2 - 10-"“) 

0.48(1.4- 10-") 

T >1 

0.84(3.5- 10-”) 

0.69 (2.8 - 10-'“) 

0.78 (2.0 - 10-") 

0.76(5.0- 10-*“) 

0.74(2.4- 10-"") 

0.74 (2.4 - 10-"") 

0.66(4.2- 10-'“) 

0.69(1.2- 10-"") 

0.50(1.6- 10-'“) 


Table 4. Spearman’s rank correlation coeficient and significance (in parenthesis) for the relation between the various bar strength proxies, with 


and without halo correction, in different morphological class bins. 


slightly larger (~ 7.5%). The bar ellipticity is fairly easy to ob¬ 
tain from photometry, therefore these equations provide an em¬ 
pirical proxy to the bar torque parameter. 

As discussed in |Martin| ( |1995[ ), e is insensitive to uniform 
changes in the relative mass of a bar of fixed shape and size. This 
constitutes one of the weaknesses of the bar strength proxies that 


exclusively rely on the bar shape (e.g. ellipticity and boxyness), 
and it partly explains the shape of the Qb - £ relation and its 
dependence on the bar-to-total ratio (as illustrated in LS2002 for 
analytical Ferrers bar models). Both Qb and account for the 
relative mass of the bar and for other stellar structures. 
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Fig. 18. Same as in Fig.[T^ but using the halo-corrected values of Qb and also the forces evaluated at the bar end (fbar). The same colour palette 
and symbols are used as in the previous plots. The Spearman’s rank correlation coefficient and significance of the relation for different T bins are 
displayed in Table 


When comparing with (right plot of Fig. 171, a bi¬ 
modality is evident. A similar trend was present in the measure¬ 
ments of |Laurikainen et ar| ( |2004a) , who separated galaxies into 
early- and late-type disks. The lower part of the plot is mostly oc¬ 
cupied by SO and early-type spirals, and the upper arm is formed 
by late-type systems. Sb and Sc galaxies are more or less evenly 
distributed in the two branches of this relation. As we argue in 
Sect. |6.3| this trend is probably mostly determined by the influ¬ 
ence of the degree of galaxy central mass concentration in the 
bar force calculation, which increases this effect more strongly 
on gb than on A™“. 


Finally, we show a monotonic rising trend between A™“ and 
E (central panel of Fig. [T7| i that resembles the Qb- e relation in 
shape, but with a convex curve. As for the bar torque parameter, 
for a narrow ellipticity interval, A™“ is larger in this branch for 
earlier types. Because the m - 2 Fourier amplitude is a good 
proxy of the bar-to-total mass ratio, this trend could again be 
partially explained by the insensitivity of e to the relative mass 
of the bar of a certain shape. 


In Fig. 18 we study the Qb ■ 


eandA““- 


Qb relations taking 


into account the halo correction and also considering an evalua¬ 
tion of the forces at the end of the bar. The observed dichotomy 
between late-type and early-type systems in these relations is 
maintained regardless of whether we use values or not. 

The relation between e and Qt evaluated at the same radius re¬ 
sembles the same trend as obtained for Qb, Likewise, the bi¬ 
modality in the comparison with A2 weakens, but is undoubt¬ 
edly present. Finally, we observe that considering the effect of 
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halos for the Qj evaluation at the bar radius obscures the segre¬ 
gation in the e - relation; it converges to a common 

trend for all Hubble types, but slightly reduces the tightness of 
the correlation. For the A2(rbar) - relation, early- 

and late-type systems still occupy different regions of the plot, 
but there is no longer an empty area separating them. 

Altogether, it appears that the bar orbits are, to a major ex¬ 
tent, controlled by the underlying disk potential, but the role of 
halos cannot be ignored for the relation between bar strength and 
shape. 


6. Characterization of bars as a function of gaiaxy 
mass 

The tendencies of the bar parameters in the Hubble sequence 
as shown in the previous section can be better understood when 
they are considered as a function of the galaxy mass. 


6.1. Bar fraction dependence on the parent galaxy mass 


As discussed in B2015, the double-peaked trend of the bar 
fraction in the Hubble sequence discussed previously might be 
linked to the bimodality found at intermediate redshifts z e 
(0.01,0.1) by|Nair & Abraham (20101. Consistent with Barazza 


|et al.| ( |2008| ), these authors found a bimodal dependence of the 
bar fraction on the central concentration and the total stellar mass 
with a minimum /bar at logio(M»/Mo) ~ 10.2. According to the 
authors, this might imply that the creation of the red and blue se- 
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Fig. 19. Bar fraction as a function of stellar mass {upper panel) and 
halo-to-stellar mass ratio {bottom panel) for the different bar detection 
criteria contemplated in this work. The error bars are binomial errors. 


quences are linked to the origin and evolution of bars, since that 
is the characteristic mass at which a bimodality in the stellar pop¬ 
ulation and fundamental properties of galaxies is observed (e.g. 


Shen et al.|2003 Kauffmann et al.|2003[|Baldry et al.|2004 l. 

Here we searched for a possible dependence between /bar 
and the total stellar mass (upper panel of Fig. [T^, but found no 
obvious trend from the visual classihcation of bars for galaxies 
of stellar masses > IO^Mq. /bar clearly drops for the faintest sys¬ 
tems ( |Sheth & S4G Team|2014| l, which is consistent with stud¬ 
ies of the bar fraction at higher redshifts (e.g. Sheth et al.|2008 
|Mendez-Abreu et al.|2010| l. 

Nevertheless, methods based on Fourier decomposition and 
ellipse fitting for the detection of bars reveal that /bar increases 
with stellar mass up to ~ 10® where it reaches a constant 

level of ~ 50% in the range lO'® - 10 "Mo with an isolated peak 
(~ 60%) at ~ IO'^'^Mq. Curiously enough, this is also the mass 
where galaxies show, statistically, the smallest amount of dark 
matter relative to the stellar mass (see again Fig. |^, and this 
mass is probably linked to the maximum /bar observed for the 


Sab galaxies in Sect. 5.1 


Cervantes Sodi et ari(2015|l found a relationship between the 


bar fraction and the halo-to-stellar mass ratio for relatively long 
bars, using the galaxy group catalogue of Yang et al. ( 2007| l and 
photometry from SDSS DR7 (stellar masses ranging from 10® ^ 
to 10*' ^ Mq): /bar drops with increasing This tendency 


Correlations: 

njar(kpc) VS. ff25..s(kpc) 

rbar(kpc) vs. /iR(kpc) 

rbar(kpc) vs. M,(Mo) 

r < 0 

0.85 (1.6- 10-“’) 

0.87(2.3 ■ 10-“’) 

0.78(1.5- 10-‘9 

T e [0,3) 

0.72(7.1 ■ 10-“) 

0 

bo 

0 

0.62(5.3- lO-'O 

T€[3,5) 

0.60(1.6- 10-") 

0.65(1.9- 10-^^) 

0.54(3.1 - 10-'') 

re [5,7] 

0.39 (5.5 ■ lO-’") 

0.53 (9.0 ■ 10-'“) 

0.12(0.0778) 

T>1 

0.68 (2.5 ■ 10-”'*) 

0.70(2.2- 10-^-*) 

0 

4^ 

0 

0 

c 


Table 5. Spearman’s rank correlation coeficient and significance for 
the relation between the bar absolute size (in kpc), the disk size (/ir and 
R 25 . 5 ) and total stellar mass. Different Hubble stage bins are considered. 


is weak when visually identifying bars (lower panel of Fig. 19 


However, it becomes clearer when relying on the Fourier and e 
lipse fitting detection of bars. This could be interpreted as caused 
by the dark matter halos stabilizing the stellar disk against bar 
formation ( |Hohl|[T97^ |Mihos et al.|1997t [DeBuhr et al.|201^ , 
which in turn would create the dependence between bar frac¬ 
tion and stellar mass (or Hubble stage) given the coupling be¬ 
tween luminous and dark matter (Fig. [^. However, we did not 
find any clear dependence of /bar on the halo-to-stellar mass ra¬ 
tio at a fixed stellar mass with any bar detection criterion (see 
Appendix]^, although this contrast could be due to our smaller 
sample size. 


6.2. Dependence of the bar parameters on galaxy mass 

It is well known that the scatter among the bar parameters is 
very high within the different Hubble type bins, and therefore 
it is natural to study the same parameters also as a function of 
galaxy mass. 

It appears that the ellipticity of a bar is independent of the 
parent galaxy mass (third panel in Fig. [20] l. The SOs are simply 
lacking the highest bar ellipticities, which might be related to 
the bar morphology in these galaxies: for example, barlenses and 
ansae are typical of early-type galaxies ( |Laurikainen et al.|201 l| l. 
For the bar torque parameter, for which the trend was obvious 
in the Hubble sequence, we observe that Qb increases towards 
lower galaxy masses. For an opposite correlation is found, 
which is expected because the trend with T was also reversed. 
None of the bar strength indexes studied as a function of the 
stellar mass is clearly segregated in the different Hubble types 
for a given galaxy mass. 

The bar sizes in physical units show a correlation with the 
total stellar mass (see the lowest panel of Fig.[20]l, which is quite 
natural (massive galaxies have larger stellar structures in gen¬ 
eral). E2005 observed that late-type spiral galaxies seem to have 
no correlation between bar size and either the disk scale length 
or the absolute magnitude of the host galaxies (with the galaxy 
parameters measured in the B band). We have reassessed this re¬ 
sult (see Table]^ and found that even though trends between rbar 
and /iR and R 2 S .5 are weaker for late-type systems, the correla¬ 
tions hold in general, regardless of T -type (but for the particular 
case of Sc and Sd galaxies we do not observe any correlation 
between rbar and M* either, as noted in E2005). Early-type spi¬ 
rals and especially SOs are the galaxies with the strongest scaling 
relations between bar and disk sizes. Eurthermore, as shown in 
the lower panels of Eig. bar-disk size relations show distinct 
trends for the systems with revised Hubble stages T > 5 on one 
hand and T < 5 on the other. 

Moreover, E2005 (see their Eig. 14) concluded that there 
is no dependence between disk-relative bar size and the parent 
galaxy absolute blue magnitude. This is more or less consistent 
with our results (see Eig. [20] l, assuming a similar behaviour be¬ 
tween blue and near-IR absolute magnitudes and limiting the 
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analysis to their T range. However, the author speculates that 
fainter galaxies should have a similar range of relative bar sizes, 
which we confirm do not. We observe a decreasing trend of the 
disk-relative bar sizes with increasing stellar masses for the late- 
type systems. On the other hand, the relative size increases with 
the stellar mass for early types. In conclusion, we find a di¬ 
chotomy in the behaviour of both the absolute and relative bar 
sizes with respect to the stellar mass of the host galaxy. 

Based on the analysis of the bar parameters, both as a func¬ 
tion of Hubble type and galaxy mass, it appears that bars in the 
SO galaxies (T < 1) are slightly different from other types. Bars 
in these galaxies are rounder and have lower Qb (which most 
probably cannot be explained solely by bulge dilution effects). 
In addition, S0“ and SO® are typically weaker based on 


Another conclusion is that the bars in the low-luminosity 
galaxies at the end of the Hubble sequence have similar ellip- 
ticities as the rest of the spirals: they do not have more oval 
shapes, for instance. In addition, while the amplitudes are 
weaker, the are higher because of the weak underlying disks 
of the parent galaxies. In the simulations by Seidel et al. ( 2015| l, 
higher values of Qb are obtained for the systems with a lower 
disk-to-total fraction, which may support the observed difference 
between Sb/c and SO/a galaxies, under the assumption that the 
former are more dominated by dark matter than the early-type 
systems (see |Bosma|1998b| for a review of the dark matter prob¬ 
lem). 


6.3. Bar strength dependence on the bulge-to-total mass 
ratio 


As the axisymmetric stellar structures are known to influence the 
measurement of bar strengths (e.g. Laurikainen et al.|2004a| l and 
the formation and evolution of bars and bulges might be con¬ 
nected (e.g. |Sheth et al.|2008[ ), we next studied in which way the 
above-derived bar strength parameters are associated with the 
properties of bulges, specifically the B/T. 


We plot in Figurej^the indicators of the prominence of the 
bar (Qb, A™™, s) as a function of B/T. The upper panel shows 
that the maximum Qb decreases with increasing B/T. We con¬ 
firm the previous result by |Laurrkainen et ar] ( |2004a) , where such 
an upper boundary was also found and interpreted as dilution of 
Qb by massive bulges. The figure also separately shows bulges 
with a Sersic index n < 2 and n > 2, as this limit is often used to 
divide them into pseudo and classical bulges, respectively ( |Ko- 
|rmendy & Kennicutt|2004t [Fisher & Drory|2008[ ). The fact that 
this upper limit exists for all kinds of bulges is consistent with 
this interpretation (only the relative mass matters). In addition, 
galaxies with bulges have a lower mean bar torque for all Hub¬ 
ble stages, as shown in Fig.[T3|(overall means of 0.23 + 0.13 and 
0.42 + 0.18 for systems with and without bulge, respectively). 
However, the discussed anti-correlation between Qb and B/T is 
stronger for bulges with n > 2 and completely disappears for 
galaxies with T > 3 (when the relative mass of the bulge, if 
present, decreases). 

The ellipticity maximum as a function of B/T shows a sim¬ 
ilar distribution as Qb. Both the upper bound and the anti¬ 
correlation only hold for bulges with high Sersic indexes, which 
can be interpreted as an effect of classical bulges making the bar 
isophotes rounder. Furthermore, a correlation is found between 
A™“ and the relative mass of the bulge if we take only galax¬ 
ies with n <2. For Sersic indexes higher than 2 the correlation 
vanishes completely. The trend is likewise strenghtend when SO 
galaxies are excluded; for these the bulge-to-total ratio and the 
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Fig. 20. Gravitational torque, bar density amplitude, bar intrinsic el¬ 
lipticity, and 1 ? 25 . 5 -relative and absolute bar size (in physical units) vs. 
stellar mass. Spearman’s rank correlation coefficient and significance 
are displayed in the plots. Different morphological types are separated 
following the colour palette shown in the legend of the first panel. The 
black symbols correspond to the running mean, and the error bars cor¬ 
respond to the standard deviation of the mean. 
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T<0 • 
T e [0-3) • 
T e [3-5] • 
T e [5-7] • 
T>7 • 


I n > 2 


p = -0.52 
p-value < 0.0001 


0.8 


1.0 



Correlations: 

rbar/R25.5 VS. Qb 

rbar/R25.5 VS. A”“ 

rbar/R25.5 VS. 6 

r < 0 

0.15 (0.2961) 

0.37 (0.0110) 

0.32 (0.0297) 

T e [0,3) 

0.36(2.3 ■ 10-") 

0.41 (2.3 ■ 10-'’) 

0.40 (2.3 • lO^") 

T 6 [3,5) 

0.55 (2.8 • 10-") 

0.53 (4.8 • 10-") 

0.43(1.9- 10-") 

T 6 [5,7] 

0.37 (6.8 • 10-') 

0.37 (7.0- 10-') 

0.17(0.0153) 

T>1 

0.09 (0.3109) 

0.18(0.0490) 

-0.13 (0.0876) 


Table 6. Spearman’s rank correlation coeficient and significance for the 
relation between the . 5 -relative bar size and the different indexes of 
the bar strength. Different Hubble stage bins are considered. 


That the Qi, value is sensitive to the bulge is one of its ad¬ 
vantages for the quantification of gravitational torques, since all 
components affect the underlying potential in the bar region and 
therefore the dynamics of the orbits making up the bar. The 


usefulness of Qi, is supported by the results reported in Seidel 
|et al.| ( |20T5] l (Fig. 12), where the photometric bar torque (as ob¬ 
tained in this work) is compared with the so-called kinematic bar 
torque (2kin), which directly measures the kinematic perturba¬ 
tion related to the bar. Qkin is derived from the radial and tangen¬ 
tial velocities extracted from the stellar velocity helds (mapped 
with the IFU spectrograph SAURON, Bacon et^|2001| l fol¬ 
lowing the recipe by [Maciejewski et al.| (|2012[). A tight corre¬ 


lation between the two independent bar strength measurements 


B/T 


Fig. 21. Bar strength indicators dependence on P4 bulge parameters. 
Spearman’s rank correlation coefficient and significance are displayed 
in each plot. The 2b, , and e as a function of the P4 bulge-to-total 

ratio are plotted in rows 1-3. Colour indicates the galaxy morphological 
type. Open symbols indicate bulges with a Sersic index n > 2, while 
filled circles correspond to n < 2. 


Sersic index in P4 might be overestimated (e.g. not all disk com¬ 
ponents such as lenses are included). 

The reported trends of the different bar strength indexes as 
a function of stellar mass can be partially explained by the de¬ 
gree of influence of the bulge on the measurements, since the 
more prominent bulges tend to reside in more massive galaxies 
(idem for other inner stellar structures that contribute to the over¬ 
all radial force held): (Qb) increases (decreases) with mass 

and B/T. With the same reasoning, B/T also determines the ten¬ 
dencies of and Qo in the Hubble sequence. In addition, the 
observed behaviour when comparing bar force measurements to 
each other (Fig. [TT]) might also be explained by bulges. Clearly, 
the effect of the bulge dilution is more pronounced in Qb, which 
can explain the T -dependent bimodality between A™“ and Qb 
and also the Qb- e relation. 


is found and is also supported by N-body simulations. In Seidel 
|et al.| ( |2015| l a sub-sample mainly composed of SOs and early- 
type S^G galaxies was used, which implies relatively prominent 
bulges and demonstrates their role in the dynamics of stars in the 
bar region. 

The formation of disky bulges in early-type galaxies can be 
driven by the evolution of bars (e.g. Laurikainen et al. 2007) 1, 
but scenarios in which pseudo-bulges form without any bar- 
dependent mechanism (e.g. stripped spirals) are also possible 
( |Laurikainen et al.|2010| l. The observed bar weakness (as mea¬ 
sured from Qb and e) that accompanies prominent bulges is most 
likely explained by the bulge dilution. Another explanation could 
be the exchange of angular momentum with gas that the bar fun- 
neled inwards and eventually formed or increased the pseudo¬ 
bulge, causing its present-day weakness. However, the correla¬ 
tions between Qb or e and B/T are weaker for low Sersic in¬ 
dexes, which does not support this interpretation. On the other 
hand, bulges formed through mergers or by the coalescence of 
giant clumps at high redshift (higher Sersix indexes) could have 
prevented bars from growing, and hence bars in the presence of 
relatively massive bulges barely evolved since they formed and 
remain weak at z = 0. A different perspective is given by consid¬ 
ering A™“ vs. B/T: the positive correlation between these two 
values seems to point in the direction of bars evolving in paral¬ 
lel with disky bulges. This link between bars and bulges is not 
observed among SOs, where a more thorough analysis of the re¬ 
lationship between bars and bulges in the Hubble sequence can 
be achieved by considering the presence of barlenses when per¬ 
forming the 2D decomposition of the surface brightness profiles 
(|Laurikainen et al.jin preparation!. 


7. Growth of bars 

A fundamental research question is whether galactic stellar bars 
become stronger, longer, thinner, and slow down as they evolve 
while exchanging angular momentum with the disk and the halo 
components (e.g. |Kormendy T979] Tremaine & Weinberg 1984| l. 
This evolution is expected from the simulation models of various 
independent groups (e.g. |Little & Carlberg||1991[ |Debattista &| 
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rA2/R25 5 



l’bar/R25.5 



l’bar/R25.5 


Fig. 22. Relation between the bar length (in kpc and scaled to the disk 
size) and the different bar strength indexes. Spearman’s rank correlation 
coefficient and significance are indicated in each plot. Colour selection 
depends on the galaxy Huhble stage. The straight line displays the linear 
fit of the cloud of points. 


Sellwoodj 19981 lAthanassoula & Misiriotis|2002[ |Athanassoula 

20031 lO’Neill & Dubinski||2003t IDebattista et al.l|2006t IValen- 


zuela & Klypin 2003 Martinez-Valpuesta et al.|2()06 1 . Observa 
tional evidence of this evolution has been investigated by differ¬ 
ent groups (e.g. lEltnegreen et al.|2007HLaurikainen et al.|2007 
|Gadotti|201 1) . 


Elmegreen et al. (2007 1 reported a correlation between ta, 
(normalized to R 25 ) and Aj , present for early and intermediate- 
type spiral galaxies based on the //-band and /fs-band mea¬ 
surements by [Laurikainen et al. ( 2006| l and Laurikainen et al. 
( |2004b| ). We also observe a definite trend between the bar size 
(both in kpc and relative to the disk size) and A™“, which gives 
evidence that longer bars are also stronger (see Eig.p^. A sim¬ 
ilar statistical trend persists for Q^, or e, but the correlation is 
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dRV3.6^m(0)[km/s/kpc] 

Fig. 23. Inner slope of the stellar component of the rotation curve as 
a function of Q^, and the normalized har length. The Spearman’s 
rank correlation coefficient and significance are displayed together with 
the linear fit of the cloud of points. For the first two panels, galaxies with 
and without massive bulges are studied separately, while for the lower 
plot colours are chosen based on T, indicating the moving average with 
black symbols. 


weaker (consistent with Laurikainen et al. 20021. We proceed by 
studying the relation between the bar sizes normalized to /? 25.5 
and the different bar strength indexes in bins of T (see the re¬ 
sults of the statistical tests in Table [^. We observe that early- 
and intermediate-type spirals, T e [0,5), are indeed the galax¬ 
ies with the tightest monotonic correlation between bar size and 
strength, independently of the bar index used. A moderate trend 
is also found among late-type spirals, which trend disappears for 
Magellanic and irregulars. There is also a moderate positive cor¬ 
relation for SOs when considering e or A™“ that practically van¬ 
ishes when taking Q^, as the bar strength proxy (for which bulges 
and barl enses may substantially change the bar strength). 


|Elmegreen et al.| ( [2007[ ) also claimed that both the bar length 
and the relative amplitude of m=2 Eourier component correlate 
with the central density, which can be understood in terms of 
galaxy evolution because galaxies with higher central rotation 
rates evolve more quickly. In Eig.j^we confirm a similar trend 
by plotting A™“ against the inner slope of the stellar contribu- 
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tion to the rotation curves, which is related to the correlation be¬ 
tween A™“ and BIT. On the other hand, high inner slope values 
are associated with low Qb, which can be understood because 
the small bar torques parameters are typical for galaxies with 
prominent stellar bulges. The anti-correlation between Qb and 
dRV3.6;im(0) appears even for galaxies without bulges. Finally, 
we reassessed the possible dependence between dRV 3 . 6 ;/m( 0 ) and 
the bar size normalized to the scale length of the disk and found 
a weak positive relationship. This trend is strengthened when the 
very late-type systems (i.e. T >7) are excluded. 

Intermediate-type spiral galaxies (e.g. She galaxies) are 
known to be the most favourable candidates to evolve secularly 
(presence of a bar, oval disk or a spiral structure; they have small 
bulges and the disk spreading is energetically favourable); they 
strip off their gas and move in the Hubble sequence towards ear¬ 
lier types ( |Kormend^|20 1 3[ ). The bar relative sizes and density 
amplitudes increase from galaxies with Sc morphology {T - 5) 
towards Sa and Sa/0 types (as shown in Figs. 11 and 12 and in 
Table [^; this is consistent with this scenario. The tightest cor¬ 
relation between the relative sizes of bars and their strengths, 
which can be understood as possible evidence of the growth of 
bars over a Hubble time, occurs for galaxies of morphological 
types Sa/O-Sc. 

For the angular momentum transfer within disk galaxies to 
be efficient, the mass of the halo near resonances is a crucial 
factor ( Athanassoula et ^|2013 l. In addition, the exchange of 
angular momentum can be hampered by dynamically hot disks 
and halos. The influence of the halo triaxiality and the dark mat¬ 
ter profile is beyond the scope of our observational analysis for 
obvious limitations. We have searched for a possible connec¬ 
tion between Mh/M» with Qb and A™™ (both depend on stel¬ 
lar mass), but found no correlation, or if one is present, it is 
very weak (Spearman’s correlation coefficient and significance 
of p = 0.18, p = 5 ■ 10“^ and p = -0.09, p = 0.048, respec¬ 
tively). However, the effect of cold gas on the evolution of bars 
is not considered in this study. The role of halos is expected to 
also be important for the evolution of bars among intermediate- 
and early-type spiral galaxies. 

The discrepant statistical behaviour of A™“ and Qb as a func¬ 
tion of r, M, , or BIT does not mean that one index is better 
than the other, as they are rather complementary. While A 2 ™” 
accounts for the relative mass of the bar, Qb directly measures 
the bar-induced tangential forces. As the tangential force de¬ 
pends on the underlying gravitation potential, it is also sensitive 
to the presence of massive spheroidals (or other vertically thick 
inner structures such as the boxy/peanut bulges) that dilute the 
Ft I < Fr > force ratio. The ellipticity of the bar, s, is also 
determined by a response of the stellar orbits to the underlying 
gravitational potential, and it therefore is a proxy for Qb- How¬ 
ever, the maximum ellipticity also depends on the morphology of 
the bar and therefore does not directly measure the bar-induced 
tangential forces. 

8. Summary and conclusions 

We analysed the bar properties of a sample of 860 face-on and 
moderately i nclined barred gal axies (i < 65°) drawn from the 
S^G survey ( Sheth et al.|||2010 i, which covers a wide range of 
galaxy masses in the local Universe. We have performed a sta¬ 
tistical analysis of the lengths and strengths of stellar bars in the 
Hubble sequence and also studied the bar fractions. 

For a sub-sample of ~ 600 barred galaxies and using the 
3.6 pm images, we have inferred the gravitational potentials with 
the polar method ( |Salo et ari|1999[ [Laurikainen & SaTo||2002t 


|Salo et al.|2010| l and calculated gravitational torques, that is, the 
ratio of the maximal tangential force to the mean axisymmetric 
radial force, to obtain an estimate of the bar strength (Qb)- In 
addition, we measured the bar strenght by the m=2 normalized 
Fourier density amplitu des (A™™) and used the bar isophotal el- 
lipticities (e) taken from |Herrera-Endoqui et al.' ( |2015[ ). 

We showed that the three proxies of the bar strength posi¬ 
tively correlate with each othe r (Fig.[T7]l, confirm ing earlier stud¬ 
ies by Block et al. P001| l and [Laurikainen et al. ( 2004a| l, among 
others. The s - Qb relation follows a curved upwards trend, 
which shows that the underlying disk potential tends to control 
the shape of the orbits making the bar. We found a segregation 
of bars determined by the Hubble type of the host galaxy. This 
might be explained by the insensitivity of e to bar mass changes, 
and this weakness of e as a bar strength index is also reflected 
in the e - A™’‘ correlation; for a fixed bar ellipticity, A™“ in¬ 
creases with T -type. In addition, we confirmed a strikingly clear 
bimodality in the relation between A™“ and Qb, which could be 
explained by the enhancement of the radial force held induced 
by axisymmetric inner structures such as bulges (bulge dilu¬ 
tion) that reduces the Qb value particularly in early-type galaxies 


(Fig.|2T]l. Indeed, both for the e-Qb and e-A™“ relations there 


seem to be two clearly distinct correlations, one associated with 
early- and intermediate-type galaxies (T < 5) and the other with 
late-type systems. We provided an estimator of the torque pa¬ 
rameter in terms of e for these two families of bars. In addition, 
we reassessed these relations by evaluating the force amplitudes 
(A 2 and Qt) at the end of the bar, finding more scattered but 
similar trends (Fig. [T8] i: this means that bulges themselves can¬ 
not explain these two T -dependent behaviours without invoking 
other physical parameters, which confirms the peculiar nature of 
late-type bars. 

We analysed the bar sizes (in kpc, as well as normalized to 
the scale length of the disk to and Ris.s) as a function of the host 
galaxy Hubble stage and total stellar mass (Figs. 11 and|20|i. We 
used bar lengths estimated visually, based on ellipse fitting, on 
the radii of maximum torque (rqb), and on the radii of maximum 
m=2 normalized Fourier amplitude (rA 2 )- All bar size estimates, 
except for that based on ellipse fitting, tend to underestimate the 
visual measurements. As in previous studies, we observed that 
bars in early-type spirals (SO/a-Sa) are typically longer than in 
the intermediate-type spirals (~Sc), and there is a decline in the 
bar size from late-type SOs towards S0“ and S0° types. How¬ 
ever, the contrast between early- and late-types is not as high as 
reported before (e.g. Martin} 19^ Erwin|2005 1. 

The richness of the sample in low-mass galaxies has al¬ 
lowed us to probe the bar properties of faint and blue systems 
at z = 0 on an unprecedented level of detail and sample size, 
even though these bars in the late-type galaxies are fairly com¬ 
plicated to measure. This is the range where we found one of the 
most conspicuous results: statistically, bars in galaxies of Hubble 
types > 5 start to become larger in size relative to the disk. This 
result is found with all of the four independent measurements 
of the bar sizes (even with rgb, which tends to underestimate 
the real bar length in late-type galaxies, as rqb moves inwards 
in galaxies with low central concentration), and regardless of the 
disk size normalization that we used. We also showed that a sim¬ 
ilar increase in the relative bar size occurs for galaxy masses 
lower than ~ IO'^’A/q. 

We checked that our quantitative estimates of bar parameters 
are consistent with the families as classified by jButa et al.| ( |20l5] l. 
We found that SB, SAB, SAB, and SAB families differ primarily 
in bar strength for all morphological types (Eigs. [TS] and [Thll. We 
also showed that strong bars are longer than weak bars. 
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According to the models of |Athanassoula| ( |2003| ) and 
Martinez-Valpuesta et al. p006|l, bars lose their angular momen¬ 


tum and become thinner and longer in time by trapping particles 
from the outer disk, which elongates their orbits or decreases 
their rotational velocity. In this scenario, the angular momentum 
left at the inner Lindblad resonance is absorbed by halo parti¬ 
cles in the surroundings of corotation and the outer Lindblad 
resonance. We found observational evidence for the evolution 
of bars, as longer bars (relative to the disk size) are also stronger 
for all the indexes (Figs. [2^. This correlation is particularly tight 
for A™“, which is the bar strength index used in the above sim¬ 
ulations. 


Furthermore, our observations showed that bars have in¬ 
creasing amplitudes and lengths from T = 5 to T = 0 
(Fig.[T§, which is consistent with the expectations of the above 
simulation models in a scenario in which late-type galaxies 
evolve secularly towards earlier types. In the same Hubble-type 
range Qb decreases towards the early-type systems, which can be 
explained by the more massive central concentrations (bulges) 
in these galaxies, which dilute the bar-induced tangential forces. 
The ellipticity of the bar is maintained nearly constant, which is 
due to the insensitivity of s to changes in the relative bar mass 
and the effect of bulges and barlenses that cause the isophotes to 
become rounder (|Martin|1995[ [Laurikainen & Salo|2002[ |Lau-| 


rikainen et al. 2004a 2007| l. All of our bar strength parame¬ 
ters decrease simultaneously only among the Hubble stages of 
r < -1 , most probably indicating that bars in these galaxies are 
less prominent than in the later Hubble stages. 


We probed possible sources of uncertainty in our Q\, mea¬ 
surements (Figs. USB an d|5): 

(1) The influence of the spiral arm torques on the Qb value. 

(2) The assumption of a constant scale height. 

(3) The effect of non-stellar emission in the bar region. 

(4) The dilution of the bar forces by the dark matter halo. 


The uncertainties associated with the first three aspects (~ 
10 - 15%) are similar in size as the estimated uncertainty re¬ 
lated to the thickness of the disk, although the mean Q\, of Sb 
and Sc galaxies can be slightly more affected by the effect of 
the spiral arms. The effect of the dark matter halo is highest for 
the latest-type systems (> 20%). The dark matter halos can re¬ 
duce the Qr value at the end of the bar by as much as 30 - 50% 
(Fig. 14 1 . Moreover, when using ICA mass maps for the Fourier 
dissection of bars, no signihcant deviation is observed in the cal¬ 
culation of the torque parameter (Fig. None of the sources 
of uncertainty alter the statistical trend of Q\, in the Hubble se¬ 
quence (Fig. [T3 ] i, although the raw bar torque parameter of in¬ 
dividual galaxies as derived in this work may deviate from the 
intrinsic bar strength for specific galaxies (e.g. those with very 
strong boxy/peanuts or strong spiral arms, massive cuspy halos, 
or high star formation rates along the bar dust lanes). 


The tendencies of the different proxies of the bar strength 
in the Hubble sequence and the impact of axisymmetric stel¬ 
lar structures on the measurements are better understood when 
studied in terms of the stellar mass of the host galaxies (Fig.p0|): 
Qb (A™“) is observed to decrease (increase) with M,, while e 
is nearly independent of M*. The behaviour of Qb and the nor¬ 
malized bar sizes of late-type spirals, Magellanic, and irregu¬ 
lar galaxies when studied as a function of M* or T is counter¬ 
intuitive: in spite of their strong Qb , they are not expected to 
be highly evolved bars. Their host galaxies are rich in gas, and 
some of them might be the result of galaxy interactions or merg¬ 
ers, especially the offset bars, but this is not sufficient to explain 
the whole picture. 
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In addition, we calculated the stellar contribution to the rota¬ 
tion curves and studied the bar properties as a function of its in¬ 
ner slope (Fig. [2^, which traces the central stellar mass concen¬ 
tration of our galaxies. We found a positive (negative) monotonic 
trend between the central rotation gradient and A 2 {Qb)- Inter¬ 
estingly, galaxies with higher central concentrations (which are 
probably more evolved) tend to host longer bars relative to the 
underlying disk size (late- and early-type system present again 
different behaviours in this relation). 


From the stellar components of the rotation curves and the 
H I maximum velocities available in the literature, we obtained a 
hrst-order estimate of the halo contribution to the rotation curve 
and of the halo-to-stellar mass ratio within the optical disk of 
our galaxies, Mh/Mt(< Ro-pi)- Assuming that the halo inside the 
disk radius forms roughly a constant fraction of the total halo 
mass, our ratios are consistent with studies based on abundance 
matching, weak lensing analysis in galaxy clusters, and halo oc¬ 
cupation distribution methods (Fig.[^. Galaxies of Hubble type 
r > 5 have a larger Mh/Mt(< /?opt)- No clear dependence of the 
bar parameters on Mh/M,(< /?opt) was found. 

We studied the bar fractions in the sample taking into 
account the different criteria for the identification of bars: 

(1) From visual classification. 

(2) By inspection of the P4 ellipticity radial profiles, with the 
detection criterion based on the presence of an ellipticity local 
maximum associated with the bar and the constancy of the posi¬ 
tion angle at that region. 

(3) From the identification of A™“, together with the constancy 
the m -2 phase, and the presence of a well-defined fourth quad¬ 
rant butterfly pattern in the torque map. 


The bar fraction shows a double-peaked distribution in the 
Hubble sequence, with a local minimum at T -4, regardless of 
the bar detection criteria (Fig. [T0| ). In comparison to the visual 
identification of bars, ellipse fitting and Fourier decomposition 
methods reduce the bar fraction by ~ 30 - 50% for galaxies with 
T > 5. This is natural because these bars are too faint to ap¬ 
pear as maxima in the ellipticity or Fourier amplitude profiles. 
Furthermore, these two methods reveal a clear monotonic trend 
between the bar fraction and both the total stellar mass and the 
halo-to-stellar mass ratio of the host galaxy, conhrming the re¬ 
sult by Cervantes Sodi et al. ( 2015| l (Fig. [T^ . However, this is 
not as clear from the visual classification. The fact that the bar 
fraction (estimated with criteria 2 and 3 above) decreases with 
Mh/M,(< f?opt) could be interpreted as a result of dark matter ha¬ 
los stabilizing the stellar disk against bar formation ( |Hohl|197^ 
|Mihos et al.|1997[ [DeBuhr et al.|201^ and might cause the de¬ 
pendence between bar fraction and stellar mass (or Hubble stage) 
given the coupling between luminous and dark matter. However, 
in contrast to |Cervantes Sodi et al.| ( |20T5] l, we did not find a sta¬ 
tistically significant difference in Mh/Mt(< /?opt) between barred 
and and non-barred systems when limiting our search to individ¬ 
ual mass bins (Fig.|E.l|i. 


Numerous N-body simulations suggest a scenario for the 
evolution of disk galaxies in which bars are expected to slow 
down as they evolve secularly, which is difficult to reconcile 
with the fact that most observed bars are fast (e.g. [Elmegre^ 
|et al.|19^|Corsini|2011| l. Recent work using direct Tremaine- 
Weingberg measurements also indicates fast bars ( [Aguerri et al.| 
|2015| l. In addition, [Rautiainen et al.| ( |2008| l showed that bars in 
early-type galaxies are typically fast rotators, while later-type 
galaxies present both fast and slow bars. In future papers we will 
reassess this topic by studying the bar pattern speed from the 
potential derived in this work, using the 2D sticky particle sim- 
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ulation method of |Salo et al.| ( |1999) t and linking the analysis to 
the bar properties obtained here. 
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Appendix A: Summary of parameters 

The following data are listed in Table [AT| (the complete version can be found at www.oulu.fi/astronomy/S4G_BARFORCE/), for 
all the 1345 disk galaxies with inclinations lower than 65° in P4 (ellipticals are excluded): 


- Column 1: Galaxy identification. 

- Column 2: Stellar mass from P3 (M»). 

- Column 3: Distance to the galaxy from P3 (D). 

- Column 4: Disk scale height (h^) estimated from the P4 disk scale length using|de Grijs|(|1998|l’s relation (or from 2MASS rk^o, 
[Speltincx et al.|2008| when /zr values were not available). 

- Column 5: Optical radius (Ropt) calculated using the P4 scale length (Ropt = 3.2/zr), or assuming Ropt = R 25.5 when /ir values 
were not available. 

- Column 6: Hi maximum velocities (yj![‘“) taken from the literature (Cosmic Flows or HyperLEDA), corrected for inclination. 

- Column 7: Maximum of the stellar contribution to the rotation curve (E^^nj)- 

- Column 8: Radius of maximum velocity of the stellar component of the rotation curve (r^g^j^^). 

- Column 9: The stellar contribution to the rotation curve evaluated at the optical radius (V 3 . 6 ;im(R opt))- 

- Column 10: Inner slope of the stellar component of the rotation curve (dRV3.6;im(0)). 

- Column 11: Halo core radius (a) computed with Eq. 15 

- Column 12: Halo velocity amplitude (Too) computed with Eq.[^ corrected to match the Hi maximum velocity. 

- Column 13: Stellar-to-halo mass ratio within the optical radius^h/4^*(< Ropt)) estimated with Eq.|l3| 


Table A.l. Stellar masses, distances, scale-heights, optical radii, maximum observed velocities, parameters of the stellar and halo components of 
the rotation curve, and halo-to-stellar mass ratios (i < 65°). 
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For the barred galaxies with measurements of Qb and rgb, the following data are listed in Table A.2 (the complete version can 
be found at www.oulu.fi/astronomy/S4G_BARFORCE/): 


- Column 1: Galaxy identification. 

- Column 2: Bar maximum gravitational torque {Qh)- Errors in Qb are associated to the highest and lowest /tR//tz values in the 
relation from |de Grijs| ( |199^ . 

- Column 3: Bar maximum gravitational torque radius (rgb). 

- Column 4: Bar maximum gravitational torque correcting for the dilution of the dark halo 

- Column 5: Bar maximum gravitational torque radius correcting for the dilution of the dark halo 

- Column 6: Bar maximum gravitational torque subtracting the effect of the spiral arms (Qb” 

- Column 7: Bar maximum gravitational torque radius subtracting the effect of the spiral arms 

- Column 8: Bar maximum gravitational torque subtracting the effect of the spiral arms and correcting for the effect of the dark 

halo ((gb”-™'y)halo-coiT^ 

- Column 9: Bar maximum gravitational torque radius subtracting the effect of the spiral arms and correcting for the effect of the 
dark halo ((r|^“-°"‘5')haio-con| 

- Column 10: Bar maximum gravitational torque using the relation from ( |Speltincx et al.|[2008| l to calculate the disk thickness 
(2b*^*)- Errors computed assigning a factor 2 thicker or thinner disk for the force calculation. 

- Column 11: Radial force profiles evaluated at the end of the bar (2T(^bar))- 

- Column 12: Halo-corrected radial force profiles evaluated at the end of the bar 


Table A.2. Gravitational torque parameters and radii, with and without spiral arms and halo correction (i < 65°). 
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For the barred galaxies with measured and (the maximum of the higher order even Fourier component amplitudes 
could only be identified for a fraction of them), the following data are listed in Table A.3 (the complete version is provided at 
WWW. oulu. fi/astronomy/S4G_B ARFORCE/): 


- Column 1: Galaxy identification. 

- Column 2; Bar intrinsic ellipticity (e). 

- Column 3; Bar maximum normalized m=2 Fourier amplitude (A™“). 

- Column 4: Bar maximum normalized m=2 Fourier amplitude radius (rA 2 )- 

- Column 5: Bar maximum normalized m=4 Fourier amplitude (A™“). 

- Column 6: Bar maximum normalized m=4 Fourier amplitude radius (rA 4 ). 

- Column 7: Bar maximum normalized m=6 Fourier amplitude (A™“). 

- Column 8: Bar maximum normalized m=6 Fourier amplitude radius (rAe)- 

- Column 9: Bar maximum normalized m=8 Fourier amplitude (A™“). 

- Column 10; Bar maximum normalized m=8 Fourier amplitude radius (rAs)- 

- Column 11: Normalized m=2 Fourier amplitude evaluated at the end of the bar (A 2 (rbar)). 


Table A.3. Bar intrinsic ellipticities and maximum normalized Fourier amplitudes and radii (m=2,4,6,8) (i < 65°). 


Galaxy 

(1) 

e 

(2) 

^max 

(3) 

rxi 

(") 

(4) 

A max 
^4 

(5) 

rA4 

(") 

(6) 

A max 
^6 

(7) 

rA6 

(") 

(8) 

A max 
^8 

(9) 

rA8 

(") 

(10) 

A2(rbar) 

(11) 

ES0012-010 

0.69 

0.48 

15.8 

0.19 

15.8 





0.32 

ES0013-016 

0.61 

0.39 

9.8 

0.13 

8.2 

0.10 

9.8 



0.26 

ESO 026-001 

0.60 

0.64 

6.8 

0.28 

6.8 

0.13 

6.8 

0.08 

8.2 

0.24 

ESO 027-001 

0.55 

0.49 

14.2 







0.48 

ESO 079-007 

0.56 

0.19 

6.8 

0.11 

8.2 





0.09 
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0.59 

0.47 

15.8 

0.25 

15.8 





0.39 

ESO 234-043 

0.66 

0.50 

14.2 







0.42 

ESO 236-039 

0.34 

0.14 

5.2 







0.10 

ESO 285-048 

0.59 

0.32 

8.2 

0.09 

12.8 





0.27 
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0.43 

0.29 

9.8 







0.22 

ESO 289-026 

0.78 

0.78 

17.2 

0.40 

17.2 
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0.19 
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0.53 
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0.33 
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0.48 

0.62 
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0.44 

ESO 345-046 
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Appendix B: Catalogue of measurements 

Figure EH summarizes all the measurements carried out in this work for the barred galaxy NGC 4314. The following plots are 
presented (for all the galaxies in our sample, an electronic version of these figures and the data used to create them can be found at 
WWW. oulu. fi/astronomy/S4G_B ARFORCE/): 

- Panel (A): 3.6 pm image de-projected to the disk plane (in units of AB mag with an upper and lower magnitude threshold of 
25 and 18). Axes are in units of arcsec. The blue line traces the visual estimate of the length and position angle of the bar. The 
light green ellipse corresponds to the isophote fitting the bar (maximum ellipticity). The purple and light blue circles (dashed 
lines) have radii equal to rA 2 and rgb, respectively. The red circumference delimits an area of radius twice the size of the bar; 

Aircie = 2- max(rvis,re). 

- Panel (B): Fourier-smoothed density with the axisymmetric component (m=0) subtracted. Contours of equal Qj are overplotted 
in white. The dotted lines indicate the regions where the tangential forces change sign, while the red line represents the visual 
measurement of the bar length and PA. The outer circle delimits a region of roughly twice the size of the bar (the same as in 
panel A). 

- Panel (C): Gravitational torque map. In the bottom of the panel the bar shows the colour scale used for the plot. Axes are in 
arcsec. Bar length and ellipticity are traced with a black solid line. The map is zoomed over the region delimited in panel (A). 
The inner dotted circle corresponds to the P4 effective radii of the bulge, in case it is fitted in P4. 

- Panel (D): Qj radial prohles for three different disk thickness models (the nominal value is the middle one), with the visual bar 
length highlighted with a dashed line. The dotted line indicates rqb. The radius range is from zero to twice the bar size. 

- Panel (E): Normalized Eourier amplitude A 2 radial profile. The dotted line indicates rA 2 - The vertical dashed line corresponds to 
the bar length. 

- Panel (E): m=2 Eourier intensity phase. 

- Panel (G): Stellar contribution to the rotation curve for three different disk thickness models. The radius is given in arcsec. The 
vertical dashed line corresponds to the bar length. The horizontal black line indicates the maximum rotation as deduced from the 
stellar mass. The solid orange line corresponds to a straight line with the slope inferred from the linear term of the m - 3 degree 
polynomial fit of the inner part of the rotation curve (dashed line). The dotted vertical orange line delimits the region where this 
fit is made, which is taken between the galaxy centre and the radius of the maximum rotation within one fourth of /? 25 . 5 - 

- Panel (FIT Summary of measurements for the galaxy with the same notation as in Tables [AT||A.2| and |A.3| 

- Panel (In Rotation curve decomposition based on the stellar mass mass rotation curve, inferred from the 3.6 pm photometry. 
The amplitude of the URC halo component (isothermal sphere) has been corrected so that the total model matches the Hi 
maximum velocity at the optical radius (horizontal green line). The vertical dashed line indicates the length of the bar. 

- Panel (J): Raw 3.6 pm and bar-only Qt profiles, with and without halo correction. The maxima of the force profiles in the bar 
region and the bar size are highlighted with vertical lines. 


Panels (I) and (J) are available only for the galaxies with usable H i velocity measurements in the literature. When the halo component of the 
RC is not estimated, the bar-only Qt profiles are included in panel (D). 
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Fig. B.l. Summary of measurements over the barred galaxy NGC 4314. 
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Fig. C.l. For a sample of galaxies, we compare the bar maximum 
normalized m = 2 Fourier amplitude {upper panel) as measured over 
3.6 ^m and ICA-corrected maps. Likewise, stellar velocities and inner 
slopes obtained from 3.6 pm {MIL = 0.53 from E2012) and mass maps 
{M/L = 0.6 from |Meidt et al.|2014| l are compared {central and lower 
panels). Error bars are computed based on boundary disk thicknesses, 
but the uncertainties coming from the MjL are indeed larger. A lower 
mass-to-light ratio at 3.6 pm, MjL = 0.47 ( [McGaugh & Schombert| 
|2014[[^uerejeta et al.|2015) l, is also tested for the comparison of max- 
imum velocities of the stellar component of the rotation curve (green 
empty circles in the central panel). The result from the statistical tests 
and the number of galaxies in the different subsamples are indicated in 
Table |C.l1 When the uncertainty in dn 1^3 6^n,(0) associated with the disk 
thickness is larger than 75 km s“' kpc“* (upper limit arbitrarily chosen), 
we exclude the galaxies from the statistics presented in this section. 


Appendix C: Measurements over P5-mass maps 

In Sect. |3.3| we quantified the effect of non-stellar contaminants 
on the bar force, finding no systematic deviation when compar¬ 
ing 2b values computed from P5-mass maps to those obtained 
directly from the 3.6 pn\ images. Using the same sample of 
34 galaxies with available ICA analysis, we repeated a similar 
test and studied whether non-stellar emission in the bar region 
can significantly change the Fourier amplitudes of bars (A 2 and 
m - 2 phase). The improvement in the obtained A™“ values 
after correcting for the non-stellar contaminants is also small 
(Fig.|^. As with Qr, noisy A 2 profiles in 3.6 pm do not smooth 
out when using P5-mass maps. 

Furthermore, we extended the study to the stellar component 
of the rotation curves, deriving the maximum velocity and the in- 
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Parameter deviation 

Mean 

Median 

|A| 

-^bin 


0.93 

0.89 

0.03 

34 

^max ^^max(ICA) 

1.04 

1.05 

0.02 

34 


1.10 

1.10 

9.64 

71 

V3.6;jin(0)/t/R ViCa(O) 

1.16 

1.04 

11.46 

51 


Table C.l. Comparison of bar strength and circular velocity measure¬ 
ments as obtained from 3.6 pm images and from ICA corrected mass 
maps for a subsample of S^'G galaxies. The mean and the median of the 
ratio of the two values are provided, together with the mean of their ab¬ 
solute difference (|A|, in the same units as in Figure [cT| l and the number 
of galaxies in the bins. 


ner slope (when possible) from P5 mass maps for a sub-sample 
of ~ 70 galaxies. In Fig. C. 1 we compare these to the original pa¬ 
rameters obtained from raw 3.6 pm images (see also Table |C.r| . 
We find that the stellar contributions to the circular velocities are 
typically overestimated by ~ 10% when the latter are used. This 
error is of the same order as the uncertainty in arising 

from the disk thickness determination. In principle, such differ¬ 
ence might be expected to arise for spirals, given the typically 
higher star formation rates (thus more emission from hot dust) at 
the spiral arm(s) region. Curiously enough, this systematic over¬ 
estimation does not seem to be strongly dependent of the stellar 
mass: it might naively be expected that the deviation is higher 
among the faintest systems (where the specific SFR is higher). In 
the case of the inner slope of Va.e/jm, where measurement uncer¬ 
tainties are larger, the mean difference is of the order of ~ 15%. 


[McGaugh & Schombert| ( |2014| l recently pointed out a lower 
mass-to-light ratio in 3.6 pm — 0.47 (with a scatter of 0.1 

dex), which is also compatible with Querejeta et al. ( |2015 1 and 
~ 15% lower than the MjL assumed in this work. In the central 
panel of Fig. jC.lj we also reassess the comparison between the 
maximum velocity of the stellar component of the rotation curve 
derived from contaminant-free maps and from 3.6 pm after as¬ 
signing f-i.Ciiim — 0.47 to the near-IR images, finding a lower 
mean difference of ~ 3%. 


In Fig. |C.2| we present the circular velocity, normalized m - 
2 Fourier amplitude, and m-2 Fourier phase obtained from the 
3.6 pm and ICA-corrected potentials for a subset of three barred 
galaxies. The force profiles of NGC 1566 and NGC0150 were 
shown in Sect. 3.3 For these galaxies, none of the measure¬ 


ments vary substantially when using P5 products. The Fourier 
profiles and circular velocity curve of NGC 4639 are also shown 
in the same figure. We have chosen this galaxy in the interest 
of presenting a case with a noticeable difference in the (stellar) 
circular velocity curve depending on whether one takes 3.6 pm 
or the ICA-corrected mass map for the gravitational potential 
calculation (while the Fourier decomposition remains very simi¬ 
lar). Still, the difference in the stellar contribution to the rotation 
curves between the mass maps and the raw 3.6 pm images is 
of the same order as the uncertainty on Vs.e/jm that is due to the 
M/L ratio. On the other hand, we have assumed a constant M/L 
for all Hubble types. However, M/L can also be expressed as a 
function of the [3.6] - [4.5] colour ( jEskew et al.|2012[|Querejeta] 
jet al.||20T5] ). In Fig. |C.2| we also indicate the resulting rotation 
curves if such a T 3 = /([3.6] - [4.5]) were taken into ac¬ 
count, obtaining a deviation that is smaller than the error bars. 
We also checked that the trends and mean differences discussed 


in Fig. C. 1 did not change when considering a colour-dependent 
M/L. 
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Fig. C.2. Circular velocity (row 1), normalized m = 2 Fourier amplitude (row 2) and m = 2 Fourier phase (row 3), inferred from the 3.6 /im 
(blue) and ICA-corrected (red) potentials for the barred galaxies NGC 1566, NGC0150, and NGC4639 (Cols. 1-3). The blue shaded area in the 
upper panels shows the uncertainty of associated to a 30% uncertainty in T 3 which is assumed to be radially constant. The red dashed 
lines correspond to P5-mass maps error bars considering the 0.1 dex uncertainty deduced in |Meidt et al.| ( (20T^ . Circular velocities computed as a 
function of a [3.6] - [4.5] colour-dependent M/L are also shown. 


Appendix D: Universal rotation curve and the halo 
correction on the force profiles 

Based on the analysis of the rotation curves of ~ 1100 spiral 
galaxies from [Persic & Salucci| ( |1995| ) and [Mathewson et al.] 
( 1992| l, PSS concluded that the shape and amplitude of an RCs 
can be described by a global parameter of the galaxy, such as the 
luminosity, and that every RC had a univocal decomposition into 
a thin exponential disk and a halo component. Counterexamples 
to this so-called universal rotation curve (URC), which was hrst 
claimed by Rubin et ar](|1980|, can be found in|Bosma](|1998a|l 


or jVerheij^ 1 1997 1 , but in spite of them, the URC remains use¬ 
ful for statistical purposes ( Battaner & Florid^|2000| l. Here we 
check how well i) the disk velocity amplitudes predicted by the 
URC ht with those derived from S'^G images, and whether ii) the 
maximum velocity amplitude agrees with the U™ from litera¬ 


ture. 


The URC predicts the rotation velocities at any radius as a 
function of the /-band absolute magnitude ([Hendry et al. 11997]!: 


yuRc(.*,Mi) = 204 ■ 10 -“ '2(a^i+22), 


Iv^ 

' ^ URC-disk 




URC-halo 


(x,Mi), 


^URC-disk(-*’“ 


(0.72 - 0.176(Mi -H 22))1.97 x'-22 

(x2 H- 0.782)1-43 ’ 


(D.l) 


(D.2) 


'^URC-haio(^>^i) = ^"[0-28 + 0.176(M, + 22)]- 
1 H- 2.25 ■ io -0 '®("i+ 22 ) 


x2 -1-2.25 ■ 10-016W+22)’ 


(D.3) 


with X = rJRopt. We note that Eq. D.2 approximates the rotation 
curve of an exponential disk, while Eq. |D.3[ corresponds to the 
isothermal halo model (see also Eq. [T5J. 

In BLS2004 the halo was modelled based on the URC pre¬ 
diction for the linear dark-to-luminous matter ratio derived by 
PSS as a function of the galaxy luminosity in the B band (LjL^, 
with L* = 10*° "^Lq and L corrected for extinction); 


= 0.4(L/L,) 


-0.9^3 1 + 1 . 5 ^(C/C*)' 


0.4 


x2-h 1.52(L/L*)0.4- 


(D.4) 


Prom the visible bulge plus disk rotation curve (Vd) eval¬ 
uated at the optical radius, BLS2004 determined the dark 
matter halo rotation amplitude by setting Mh/M,|uRc(4c) = 
yuRC-haio( 4 i:)VVd(x)^ at X = 1 . 

Taking advantage of the availability of H i data for most of 
the galaxies in the S'^G sample, we test how well the predicted 
velocities (based on the URC) match the observed circular ro¬ 
tation (Pig. |D3. In the statistical sense, the URC model esti¬ 
mates the RC maximum velocities fairly well. We also compare 
the URC disk component with our measured Vs.e/im at their peak 
values, observing disagreement among the faintest galaxies. 
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value (deviation < 10%). Moreover, we searched for a possible 
relation between the disk inclination and Mh/M*, finding no de¬ 
pendence and thus excluding a possible bias due to observation 
uncertainties for highly inclined and face-on galaxies that might 
affect our halo corrections. We discarded any major underesti¬ 
mation of the halo amplitude (and thus of the halo correction) 
that would be due to the assumption of a radially constant mass- 
to-light ratio (e.g. Portinari & Salucci 2010[ ) because this was 
one of the factors that determined the nominal 30% uncertainty 
(or even higher, according to |Meidt et al.|2014) l on the mass-to- 
light ratio estimation by E2012, and radial colour gradients are 
not that pronounced in our wavelengths. 


w 
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CO cr 
E => 
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Fig. D.l. First row: Maximum H i circular velocity from the literature 
(corrected for inclination) vs. maximum of the universal rotation curve. 
Second row: Ratio of the maximum of the circular velocity disk compo¬ 
nent as derived from the URC and from the 3.6 pm potential vs. absolute 
corrected /-band magnitude. 


In Fig. D.2 our rotation curve decomposition models are 
compared with Ha rotation curves from Daigle et al.| ( |2006| l 
and Dicaire et al. (|2008 1 , for a sub-sample of SINGS galaxies 
( Kennicutt et al.|2003| l. Likew ise, in Fig. |D.3|we s tudy the ob¬ 
served HI velocity curves from |de Blok et al. ( 2008| l, taken from 
THINGS fWalter et al.|2008l ). Altogether, we have ten galaxies 
(three of which are barred) to compare to our ad hoc RC mod¬ 
els. In almost all the galaxies, the modelled curve seems to fit 
the observations fairly well, the outliers being NGC2841 (with 
a rotation curve of similar shape as Vs.e/jm) and the inner part of 
the rotation curve of the barred galaxy NGC2903. For all the 
systems, we obtain a ~ 15 - 20% mean dispersion of our models 
with respect to the observed rotation curves across the optical 
disk, which could be due to the gas component lacking in our 
decompositions. 

Additionally, in the lower panel of Fig. [7] we tested the ro¬ 
bustness of the halo model and the correction of the forces by 
assessing the effect of a measurement error on the used input pa¬ 
rameters. We assumed a 20% overestimation error on the mass- 
to-light ratio, the distance to the galaxy, and the observed veloc¬ 
ity; and we tested the effect of reducing by 1 mag the /-band total 
luminosity. We also checked the effect of matching the observed 
maximum velocity at 2.2/!r (radius of the maximum velocity of 
an exponential disk; |Freeman||1970|l instead of at /?opt- A pos¬ 
sible inaccuracy assuming that the maximum circular rotation 
takes place at the optical radius was also considered when testing 
the overestimation of the observed velocity (e.g. the difference 
between either the absolute maximum or the asymptotic limit 
of VuRC and its value at /?opt is always lower than 20%). None 
of these variations result in a significant change in the 
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Fig. D.2. Rotation curve decomposition models for a sub-sample of 4 galaxies from the SINGS sample (Kennicutt et al. 2003) with disk inclinations 
lower than 65° and the Ha rotation curves from Daigle et al. (2006) and Dicare et al. (2008). The RCs have been corrected for inclination using 
P4 orientations. The legend in the first pannel shows the meaning of the different lines and points. For every pannel, the x-axis indicates the 
galactocentric radius in units of arcsec (limited to the optical radius, ~ 3.2/2 r) and the y-axis corresponds to the rotation velocity in units of km/s. 
The shaded area shows the uncertainty in our rotation curve models due to a 30% uncertainty in the M/L ratio. Galaxy RCs are sorted in terms of 
their mass (increasing from left to right, top to bottom). 
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Fig. D.3. The same as in the previous figure but using the Hi rotation curves from de Blok et al. (2008) for a sub-sample of 6 galaxies taken from 
the THINGS sample (Walter et al. 2008). 
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Appendix E: Halo-to-stellar mass ratio and the bar 
fraction 


In Sect. 3.6.1 we showed a clear dependence between the bar 
fraction and both the halo-to-stellar mass ratio, evaluated inside 
the optical radius, and the stellar mass of the host galaxy. We also 
showed that Mh/M,(< Ropt) scales inversely with M, for masses 
M* < 10"Mo. To probe co-dependences between M* 

and /bar, we studied the distribution of barred and non-barred 
galaxies (as classified in |Buta et al.|2015| l as a function of Mh/M* 
for a given M, bin. We observe no significant difference be¬ 
tween them (see Fig. ED- To confirm this result, we performed 
a two-sample Kolmogorov-Smirnov (see the p-values listed in 
Table |^, finding that the cumulative distribution function of 
Mh/M* is similar for barred and non-barred galaxies, given a 
certain stellar mass range. 

We also performed K-S tests for galaxies with relative bar 
sizes smaller or greater than the median rbar/ZiR and also for the 
different bar families (SAB and SAB on one hand, and SB and 
SAB on the other), all of them compared to the family of non- 
barred galaxies. In all cases, the cumulative distribution func¬ 
tion of Mb/M» is not significantly different for barred and non- 
barred galaxies (regardless of size, strength, and bar detection 
criterion), given a certain stellar mass range. 
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Fig. E.l. Co-dependence between halo-to-stellar mass ratio, the total 
stellar mass, and the bar fraction (B2015) are studied by binning the 
MhlM,-M, space. In the upper panel, red and blue filled circles corre¬ 
spond to barred and non-barred galaxies, respectively. For a given stel¬ 
lar mass interval, we display the distribution of barred and non-barred 
galaxies vs. Mb/A/,. P-values after performing a two-sample K-S test 
are listed in Table |E. i'l l: according to the test, the differences in the dis¬ 
tribution of barred and non-barred galaxies are not statistically signifi¬ 
cant for any of the mass intervals. The green lines correspond to the total 
halo-to-stellar ratio from |Moster et al. l |2010| l scaled down by constant 
factors ~ 0.04 (solid line), ~ 0.1 and ~ 0.01 (dashed lines). 
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1 K-S p-value: 

8.5 < log,o(M./Me) < 9 

9 < log,o(M./Mo) < 9.5 

1 9.5 < logio(M./Mo) < 10 

10<logio(M*/Mo)< 10.5 1 

10.5<log,o(M./Mo)< 11 1 

Visual 

0.63 

0.17 

0.01 

0.25 

0.55 

Ellipse fitting 

0.34 

0.02 

0.32 

0.25 

0.11 

3 & !2b 

0.81 

0.03 

0.15 

0.63 

0.14 




Tbar > 0.9 /ir 



Visual 

0.36 

0.28 

0.04 

0.16 

0.29 

Ellipse fitting 

0.34 

0.20 

0.47 

0.09 

0.09 

3 & Qb 

0.68 

1.00 

0.38 

0.19 

0.18 




rbar < 0.9 /ir 



Visual 

0.52 

0.26 

0.01 

0.69 

0.46 

Ellipse fitting 

0.43 

0.03 

0.54 

0.68 

0.14 

3 A”" & gb 

0.00 

0.00 

0.08 

0.74 

0.21 




SB+SAB 



Visual 

0.85 

0.11 

0.07 

0.11 

0.27 

Ellipse fitting 

0.93 

0.02 

0.42 

0.10 

0.36 

3 & Qb 

0.99 

0.04 

0.44 

0.29 

0.44 




SAB + SAB 



Visual 

0.09 

0.47 

0.00 

0.67 

0.22 

Ellipse fitting 

0.21 

0.09 

0.31 

0.72 

0.01 

3 & 0b 

0.45 

0.35 

0.01 

0.80 

0.03 


Table E.l. Significance level of the Kolmogorov-Smirnov statistic between two sub-samples of barred and non-barred galaxies studied as a 
function of the stellar-to-halo mass ratio. The sub-samples were drawn from the stellar mass bins shown in the lower panel of Figure (the M, 
thresholds are displayed in the first row). K-S tests are also performed for galaxies with relative bar sizes lower or greater than the median 
and also for the different bar families (SAB and S AB on one side, and SB and SAB on the other), all of them compared to the family of non-barred 
galaxies. In general, the distribution of Mh/M, is not significantly different for barred and non barred galaxies (regardless of size and strength), 
given a certain stellar mass range (indicated by p-values > 0.01). 
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